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1.  Introduction 

This  research  program  focuses  on  the  resuscitation  of  traumatic  casualties  in  hypovolemic  shock. 
The  principal  objectives  are  to  determine: 

•  when  it  is  necessary  for  a  medic  on  the  battlefield  to  resuscitate  a  trauma  victim, 

•  how  the  effectiveness  of  the  resuscitation  effort  can  be  measured  objectively  with  outcome 
predictors,  and 

•  what  is  the  desired  end  point  of  resuscitation:  can  it  be  defined  empirically  by  the  survivors’  values? 
Human  studies  at  a  major  Level  I  trauma  center,  Los  Angeles  County+USC  Medical  Center 
(LAC+USCMC),  serve  to  precisely  define  the  pathophysiology  of  traumatic/hemorrhagic  shock  and  test 
the  effectiveness  of  noninvasive  (or  minimally  invasive)  sensors  of  tissue  perfusion  failure.  The 
investigation  entails  the  consolidation  of  three  large  existing  trauma  data  sets  at  LAC+USCMC  into  a 
new  database  containing  existing  medical  record  data  without  the  possibility  of  patient  identification. 
Our  use  of  the  term  “existing  medical  record,  EMR”  data  in  this  report  automatically  carries  the  qualifier 
“without  the  possibility  of  patient  identification”.  The  current  database  tracks  the  trauma  victim 
beginning  at  the  emergency  department,  then  to  the  operating  room,  the  radiology  suite  when  indicated, 
and  the  surgical  intensive  care  unit  (SICU).  It  provides  an  excellent  reference  source  and  testing  facility 
for  the  current  state-of-the-art  in  U.S.  resuscitative  medicine.  Our  goal  is  to  make  the  entire  trauma 
database  accessible  to  the  combat-casualty-care  community  in  an  easy-to-use,  straightforward  fashion. 
Thus,  each  member  of  the  community  will  be  able  to  derive  his  or  her  own  conclusions  from  a 
comprehensive,  objective,  unbiased  database.  Taken  as  a  whole,  the  new  database  will  allow  optimum 
resuscitation  strategies  to  be  developed  for  the  Army  medic  and  will  serve  as  a  guide  for  implementing 
new  sensor  strategies.  The  LAC+USCMC  data  is  being  “packaged”  with  a  database  management 
program  and  augmented  with  a  user-friendly,  menu-driven  program  for  standard  analyses  and  data 
manipulation. 

2  Body  of  Report 

During  this  first  grant  year,  particular  attention  was  paid  to  building  the  infrastructure  for  the 
trauma  database  and  populating  the  database  with  EMR  patient  case  histories.  In  addition,  outcome 
predictors  were  investigated  within  the  context  of  a  mathematical  “nearest  neighbor”  model  as  well  as  in 
studies  involving  single  parameter  predictors. 

2.1  Human  Hemodynamic  Data  Processing,  Analysis,  and  Interpretation 

Despite  the  large  number  of  existing  trauma  cases  in  the  William  C.  Shoemaker  (WCS)  database 
at  LAC+USCMC,  it  became  clear  in  the  early  stages  of  this  study  that  more  retrospective  cases  were 
needed.  In  part,  this  is  because  the  sensors  and  methods  of  noninvasive  hemodynamic  data  acquisition 
have  improved  with  time.  As  a  result,  the  more  recent  data  sets  obtained  within  the  past  five  years  are 
substantially  higher  in  quality  than  much  of  the  data  acquired  prior  to  1996.  (Approximately  1500 
trauma  cases  were  documented  at  hospitals  other  than  LAC+USCMC  prior  to  1996.)  Thus,  the  recent 
LAC+USCMC  data  allows  us  to  address  the  questions  of  Section  1  in  a  more  comprehensive  fashion. 
At  the  same  time  the  latest  data  has  proven  very  useful  in  studies  of  outcome  predictors.  Both  the 
stochastic  prediction  methodology  [Ref  1]  and  the  search  for  single  parameter  predictors  yielded 
significant  results  with  the  latest  data  sets.  The  acquisition  of  more  data  allows  the  phase  space  of  the 
mathematical  model  to  be  more  fully  populated  and  therefore  improved.  In  addition,  both  areas  of 
inquiry  greatly  benefit  from  the  existence  of  a  large  validation  data  set.  In  summary,  there  is  a  strong 
rationale  for  acquiring  more  data  from  existing  medical  records  for  inclusion  in  the  trauma  database. 
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During  the  period  September  29,  2001  through  September  28,  2002,  approximately  180  new 
trauma  cases  were  added  to  the  WCS  database.  This  is  an  increase  of  about  40%  compared  to  the 
previous  year.  Funds  from  the  current  Army  grant  were  used  for  the  processing,  analysis,  and 
interpretation  of  the  2002  data  as  well  as  all  other  retrospective  data  acquired  between  1996  and  2001. 
Methods  for  noninvasive  data  acquisition  are  described  below,  and  the  results  of  year  1  efforts  to 
analyze  and  interpret  the  WCS  database  are  discussed.  This  investigation  was  performed  outside  the 
framework  of  the  Main  Project  Database  because  the  Database  was  concurrently  constructed  while  the 
hemodynamic  research  was  being  performed.  A  single-parameter  predictor  study  was  initiated  after  the 
Database  was  completed.  This  is  discussed  in  Section  2.2.8. 

2.1.1  General  Methods  for  Noninvasive  Hemodynamic  Data  Acquisition. 

At  Los  Angeles  County  +  USC  Medical  Center  (LAC+USCMC),  electrodes  are  attached  to  the 
base  of  the  neck  and  to  the  lower  lateral  chest  wall  for  thoracic  electrical  impedance  measurements  of 
cardiac  output,  and  a  finger  probe  is  used  for  pulse  oximetry  (SapCh).  Also  monitored  are  arterial  blood 
pressure,  heart  rate,  and  ECG  as  general  circulatory  measures.  In  addition,  physicians  evaluate  tissue 
perfusion/oxygenation  by  transcutaneous  oxygen/C02  sensors  that  have  the  same  Clark  polarographic 
oxygen  sensor  used  routinely  in  the  standard  blood  gas  analyzer.  While  the  oxygen  tension  of  heated 
skin  of  the  chest  does  not  reflect  the  state  of  oxygenation  of  all  tissues,  it  has  been  shown  to  be  a 
sensitive  early  warning  of  hypovolemic  shock  [Refs.  1-3].  All  parameters  are  monitored  and  recorded 
continuously  after  emergency  department  (ED)  admission  or  following  the  onset  of  surgical  operations. 
In  so  far  as  is  possible,  LAC+USCMC  specifically  records  these  hemodynamic  variables  before  and 
after  each  therapy,  beginning  with  admission  to  the  ED,  following  the  patient  to  the  operating  room 
(OR)  and  intensive  care  unit  (ICU).  Each  therapeutic  intervention  is  noted  in  time,  dose,  duration  of 
infusion,  and  time  after  admission.  Each  therapy  is  given  one-at-a-time,  if  possible  (except  for  clinical 
exigencies),  with  hemodynamic  monitoring  before,  during  and  after  its  administration.  Perioperatively, 
noninvasive  hemodynamic  data  in  high-risk  surgical  patients  is  started  in  the  OR  or  ICU  prior  to 
anesthesia  and  maintained  throughout  the  operation.  Blood  losses  are  estimated  by  vital  signs, 
hematocrits,  urine  output,  sponge  counts,  amount  of  fluid  in  the  suction  bottles,  number  of  transfusions 
needed  to  maintain  normal  hematocrit  values,  and  clinical  observations.  Estimated  blood  losses  are 
promptly  replaced. 

2.1.2  Analysis  and  Interpretation  of  the  Database  in  Terms  of  Outcome  Prediction 

Over  the  past  year,  we  confirmed  that  multiple  noninvasive  monitoring  systems  are  feasible  and 
highly  useful  methods  to  acquire  early  information  from  patients  with  acute  life-threatening  trauma  [Ref. 
1].  These  noninvasively  monitored  data: 

•  were  sufficiently  close  to  the  information  supplied  by  the  PA  catheter  to  be  a  useful  surrogate  for 
invasive  monitoring; 

•  could  be  used  immediately  after  Emergency  Department  (ED)  admission  and  throughout  the  hospital 
(analogous  to  use  in  prehospital  settings), 

•  could  be  used  to  compare  the  time  course  of  clinical  and  hemodynamic  parameters; 

•  provided  continuous,  real  time  displays  and  data  streams  rather  than  traditional  methods  that  provide 
infrequent  and  intermittent  “snapshots”; 

•  were  utilized  to  describe  the  physiologic  and  pathologic  patterns  of  patients  with  severe  trauma  and 
other  emergencies  during  resuscitation  immediately  after  admission  to  the  ED; 

•  revealed  hemodynamic  patterns  of  severe  trauma  that  are  related  to  outcome;  and 

•  most  importantly,  provided  information  to  a  user-friendly  program  that  can  be  used  effectively  by 
nurses  and  other  paramedic  personnel  at  the  bedside  of  acutely  ill  patients. 
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In  addition  to  the  conventional  vital  signs,  hematocrit,  urine  output,  blood  gases,  and  other 
clinical  evaluations,  we  are  currently  evaluating  noninvasive  hemodynamic  monitoring  used  with  a  new 
mathematical  analysis  for  application  to  acute  high-risk  critical  illness,  shock,  and  potentially  lethal 
organ  failure  beginning  with  ED  admission.  To  evaluate  the  accuracy  of  the  random  (stochastic) 
analysis,  we  are  investigating  various  subsets  of  trauma  including  blunt  or  penetrating  injury  to  chest, 
abdomen,  extremity,  head,  neck  and  spinal  cord  in  various  combinations. 

We  plan  to  employ  hemodynamic  patterns  to  develop  outcome  predictors  using  the  new 
mathematical  analysis  and  data  mining  programs  to  elucidate  underlying  hemodynamic  deficiencies,  to 
evaluate  and  compare  the  effectiveness  of  various  therapies  in  previously  recorded  databases,  and  to 
display  possible  optimal  therapy  in  a  real-time  mode  for  each  new  individual  patient  at  each  stage  of  the 
patient’s  hospital  course.  In  addition  we  plan  to  utilize  additional  tools  to  evaluate  their  concordance 
with  established  patterns  of  monitored  data  and  physiologic  patterns.  Therapeutic  policies  based  on 
analyses  of  the  databases,  along  with  outcomes  of  other  competing  therapeutic  and  diagnostic 
modalities,  are  also  being  investigated  [Refs.  2-4]. 


As  shown  in  Table  1,  we  have  found  markedly  different  misclassification  rates  in  the  same  series 
of  trauma  patients  from  single  isolated  values  of  vital  signs,  cardiac  index,  APACHE  n,  discriminant 
analysis,  and  now  by  using  the  present  version  of  the  stochastic  analysis  program.  Heart  rate,  mean 
arterial  pressure  (MAP),  and  cardiac  index  values  were  the  initial  or  lowest  values  as  stated,  APACHE  II 
values  were  calculated  on  a  daily  basis,  discriminant  analysis  covered  the  resuscitation  period,  and  the 
probability  analysis  was  calculated  for  each  set  of  values;  the  first  set  of  values  are  presented  here.  The 
data  show  the  limitations  of  these  single  values.  By  contrast,  there  are  advantages  of  a  computerized 
database  of  multiple  noninvasive  variables  organized  over  time. 


Table  1  Misclassifications  in  Outcome  Prediction 
by  Evidence-based  Analyses 


Method 

Initial  Heart  Rate 
Initial  MAP 
Lowest  MAP 
Initial  Cardiac  index 
APACHE  II 
Discriminant  Analysis 
Probability  Analysis, 


Criteria 

S  <95,  NS>96  beat/min 
S  >85,  NS<70  mmHg 
S  >50,  NS<  50  mmHg 
S>3.8,  NS<  3.8  L/min/m2 
S  <27,  NS  >27 
(See  Ref.  22) 

S>75%,  NS  <55% 


Misclassification  Rate 
(70/159)  45% 

(76/159)  47% 
(83/159)  52% 
(72/159)  43% 

(30/97)  31% 

(23/151)  15.2% 
(16/165)  9.2% 


Initial  studies  showed  over  90%  accuracy  in  predicting  outcome  in  the  first  12  hours  after  ED 
admission  using  the  new  mathematical  (nearest  neighbor)  analysis  and  display  program  that  estimates 
the  probabilities  of  survival  based  on  the  primary  diagnosis,  selected  relevant  covariates,  and 
noninvasive  hemodynamic  monitoring  of  cardiac,  respiratory,  and  tissue  perfusion  functions.  We  plan 
to  study  approximately  150  severe  high-risk  patients  during  the  coming  year.  They  are  drawn  from  a 
pool  of  over  5000-trauma  admissions/year  to  our  Trauma/Critical  Care  Service.  We  plan  to  continue  our 
study  of  the  most  severely  ill  of  these  trauma  patients. 


The  monitoring  and  information  system  also  displays  a  decision  support  system  to  provide  the 
physician  at  the  bedside  with  a  potentially  useful  method  for  optimizing  the  therapy  by  evaluating  each 
potential  therapy  based  on  responses  of  the  nearest  neighbors  in  terms  of  reduced  mortality.  The  system 
may  be  used  to  diagnose,  evaluate,  and  suggest  therapy  for  the  early  hemodynamic  aspects  of  circulatory 
dysfunction,  shock,  and  critical  illness  anywhere  in  the  hospital  including  the  ED,  operating  room,  post 
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anesthesia  recovery  room,  intensive  care  unit  (ICU),  step-down  units,  hospital  floors,  pre-hospital  areas, 
and  outpatient  clinics.  The  major  underlying  challenge  is  to  develop  an  information  system  based  on  a 
sufficiently  large  noninvasive  database  with  adequate  numbers  of  nearest  neighbors  to  predict  outcome. 
Ultimately  this  data  base  will  support  a  therapeutic  decision  support  system  that  will  display  the  effects 
of  therapies  that  had  been  given  to  similar  patients  (nearest  neighbors),  and  identify  the  therapies  that 
carry  the  greatest  expectation  of  optimizing  survival.  Since  the  computational  burden  does  not  appear  to 
be  intense,  the  system  was  designed  to  function  in  a  real-time  mode  at  the  bedside.  Alternative  therapies 
and  the  probability  for  reduced  severity  (mortality)  will  be  displayed  at  each  stage  of  acute  life- 
threatening  illness  as  well  as  before  and  after  each  therapy.  These  data  along  with  all  supporting 
evidence  will  provide  the  physician  with  maximal  information  for  selecting  therapies  to  optimize 
outcome.  The  quest  for  early  diagnosis  and  therapy  of  circulatory  shock  is  analogous  to  that  of  cancer, 
where  early  therapy  is  far  more  effective  and  less  costly.  In  fact  the  concept  of  the  “Golden  Hour”,  the 
major  descriptor  of  the  response  to  trauma,  is  largely  based  on  the  need  to  provide  therapy  in  the  earliest 
stage  of  traumatic  injury.  In  preliminary  studies  on  a  consecutively  monitored  group  of  severely 
traumatized  patients  with  30%  mortality,  the  stochastic  analysis  program  correctly  identified  survivors 
and  nonsurvivors  with  over  90%  accuracy  in  the  first  12  to  24  hours  after  admission.  There  were  16/165 
patients  or  9.7%  who  were  misclassified  in  this  initial  resuscitation  period  [Ref.  3].  Inspection  of  the 
data  revealed  that  accuracy  might  be  limited  by  the  number  of  patients  that  comprised  the  pool  of 
nearest  neighbors  for  each  etiologic  category  in  the  database.  We  will  address  this  problem  by 
increasing  the  number  of  trauma  subsets  and  other  diagnostic  categories  studied. 

2.2  Development  of  the  Trauma  Database  With  Emphasis  on  Hemodynamics  Including  Blood 

Flow  and  Adequacy  of  Tissue  Perfusion 

As  part  of  this  project,  we  are  assembling  a  large  trauma  database  using  data  gathered  at  Los 
Angeles  County  +  USC  Medical  Center  (LAC+USCMC),  with  the  capacity  for  additional  data 
integration  from  other  sources.  The  goal  is  to  produce  a  comprehensive,  objective,  unbiased  database 
that  can  be  independently  mined  by  members  of  the  combat-casualty-care  community.  The  database  is 
expected  to  provide  insight  into  the  resuscitative  outcome  of  combat  casualties  and  civilian  trauma 
victims.  One  of  the  major  goals  of  this  project  is  to  identify  physiological  parameters  that  are  strong 
predictors  of  outcome  or  reflect  the  need  for  further  therapy. 

Our  original  statement  of  work  indicated  that  there  were  only  two  databases  that  needed  to  be 
combined  as  part  of  this  project.  These  were  the  William  C.  Shoemaker  (WCS)  database  and  the  SICU 
database.  However,  it  quickly  became  apparent  that  the  benefits  of  including  a  third  database,  the 
trauma  registry  (TR),  were  very  large.  The  TR  contains  the  official  LAC+USCMC  summary  of  each 
trauma  case  and  includes  admission  and  discharge  data.  (See  Table  8  for  a  partial  parameter  list.)  Its 
role  is  to  improve  the  efficiency  of  the  trauma  care  system  [Ref.  5].  TR  records  at  LAC+USCMC  are 
collected  into  a  relational  database  that  is  amenable  to  computer  processing.  As  a  result,  The  Main 
Project  Database  currently  consists  of  three  dissimilar  elements:  the  WCS  database,  the  SICU  database, 
and  the  TR.  The  WCS  database  tracks  patients  using  predominantly  non-in vasive  measures  beginning 
with  resuscitation  in  the  emergency  department  (ED)  through  the  RAD/OR  and  into  the  SICU.  Patient 
histories  terminate  at  discharge  or  death.  Approximately  20-30  physiological  parameters  are  recorded 
for  -689  individual  patients  from  1996  to  present.  Most  measurements  are  non-in  vasive  and  the  quality 
of  these  data  is  very  high.  The  catalog  of  WCS  parameters  has  evolved  because  of  the  introduction  of 
new  non-in  vasive  sensors.  However,  the  database  is  anchored  by  a  set  of  common  parameters.  The 
massive  SICU  database  provides  a  comprehensive  electronic  record  of  all  customary  diagnostic 
measurements,  procedures,  and  therapies.  Many  are  invasive  in  nature,  and  most  are  considered  “gold 
standard”  measurements.  Supplemental  data  is  supplied  by  the  TR,  including  information  obtained  by 
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paramedics  and  ED  personnel.  The  tools  being  developed  to  probe  the  database  are  flexible  and  are  not 
limited  to  specific  methodologies  or  techniques. 

2.2.1  Database  Assembly 

As  noted  above,  the  database  constructed  as  part  of  this  project  serves  as  a  tool  for  assessing  a 
variety  of  interventions  and  is  expected  to  provide  insight  into  the  resuscitative  outcome  of  combat 
casualties  and  civilian  trauma  victims.  One  of  the  major  goals  of  this  project  is  to  identify  physiological 
parameters  that  are  strong  predictors  of  outcome  or  reflect  the  need  for  further  or  modification  of 
therapy.  At  this  point,  we  don’t  know  whether  or  not  a  sensor  system  can  fully  replace  a  medic’s 
judgment.  It  is  most  likely  that  new  sensors  or  arrays  of  sensors  will  be  identified  that  can  enhance  the 
ability  of  the  medic  or  other  responder  to  asses  the  status  of  a  casualty.  Some  predictors  may  require  a 
combination  of  sensor  readings  and  a  medic’s  judgement.  Other  questions  related  to  the  most  desirable 
beginning  and  end  points  of  the  resuscitation  process  and  what  therapies  produce  the  most  survivors  can 
be  explored  with  the  existing  database.  Overall,  the  data  mining  process  is  quite  general  and  is  not 
restricted  to  any  single  methodology.  Many  distinct  mining  strategies  are  expected  to  evolve  with  time. 

2.2.2  The  “Toy  Database” 

During  the  current  grant  year,  a  test  data  set  referred  to  as  the  “toy  database”  was  initially 
constructed.  It  consists  of  small  portions  of  the  William  C.  Shoemaker  (WCS)  database,  the  SICU 
database,  and  the  LAC+USCMC  trauma  registry  (TR).  The  toy  database  was  assembled  to  test  the 
practicability  of  combining  the  three  diverse  databases.  In  the  case  of  the  SICU  database,  the  number  of 
extracted  parameters  was  limited  to  ~30.  Approximately  24  parameters  were  extracted  from  the  TR,  and 
all  data  from  the  WCS  Heart  Rate  Variability  (HRV)  study  (discussed  below)  were  included.  The  WCS 
database  takes  the  form  of  entries  in  Microsoft  Excel  spreadsheets.  Both  the  SICU  and  TR  databases  are 
relational  in  nature  and  are  hosted  by  two  different  database  programs,  Sybase  (SICU)  and  Foxbase 
(TR),  respectively.  The  objective  behind  assembling  the  toy  database  was  to  identify  problems  involved 
in  the  consolidation  of  these  three  data  sets.  The  ICU  database  is  by  far  the  largest;  it  is  50  to  100  times 
larger  than  the  WCS  database.  The  TR  is  in  itself  a  very  large  database,  but  the  amount  of 
data/information  that  is  extracted  for  individual  trauma  patients  is  approximately  equal  to  the 
corresponding  patient  file  sizes  in  the  WCS  database.  Data  dictionaries  had  to  be  developed  for  all  three 
data  sets  because  they  had  never  been  used  outside  of  their  targeted  application.  Figure  1  below 
illustrates  the  intersection  of  the  three  databases;  the  area  with  the  star  represents  the  data  records  that 
became  part  of  the  toy  database.  A  total  of  53  patient  histories  were  common  to  all  three  data  sets 
(WCS,  SICU,  and  the  TR).  These  were  merged  together  into  a  single  database  using  the  Sybase 
Anywhere  8.0  database  program.  A  listing  of  toy  database  parameters  is  provided  in  the  Table  2  below. 
This  new  test  database  was  probed  and  queried  primarily  with  SQL  commands  to  verify  that  it  was 
properly  ordered  and  complete. 

It  was  found  that  a  considerable  effort  is  involved  in  sorting  and  bringing  large  quantities  of 
diverse  data  from  the  SICU  database  into  the  main  Sybase  database.  In  addition,  there  is  a  strong 
preference  for  data  obtained  at  the  earliest  times  following  patient  admission  at  LAC+USCMC.  The 
SICU  data  typically  begins  6  to  24  hours  after  admission,  whereas  the  WCS  data  typically  begins  at 
admission.  As  a  result,  our  initial  efforts  focused  on  merging  the  WCS  database  and  the  TR  database 
into  a  single  project  database.  This  allowed  the  early  time  hemodynamics  of  trauma  patients  to  be 
examined  in  detail  during  the  first  grant  year.  At  the  request  of  the  Government,  all  relevant  comments 
(i.e.  character  strings)  in  the  WCS  data  sets  were  parsed  so  that  they  could  be  readily  interpreted  as  part 
of  the  Main  Project  Database.  This  was  a  significant  task  involving  the  parsing  of  more  than  50,000 
comments.  This  task  required  two  and  a  half  months  to  complete. 
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W.  C.  Shoemaker  Database 


Figure  1 .  Data  for  the  Toy  Database  was  extracted  from  the  intersection  of  three 
large  data  sets. 


Table  2.  Parameters  Used  in  the  Toy  Database 


Table:  HRV_COTD  -  HRV  Data  (WCS) 


Column  Name 

Format 

Description 

dt 

Timestamp 

Date/Time  of  measurement 

hrv_psn 

tinyint 

Pt  number  from  HRV  study 

hrv_evt 

integer 

Event  number  from  HRV  study 

cotd 

numeric(4,2) 

Cardiac  Output  (Thermodilution)  value 

Table:  HRV.FRQ-H 

RV  Data  (WCS) 

Column  Name 

Format 

Description 

dt 

Timestamp 

Date/Time  of  measurement 

hrv_psn 

tinyint 

Pt  number  from  HRV  study 

hrv_evt 

integer 

Event  number  from  HRV  study 

lfa 

double 

Low  Frequency  Area 

rfa 

double 

Respiratory  Frequency  Area 

fract 

double 

LFA/RFA 

mhr 

double 

mean  HR 
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Table:  HRVJNTEG- 

HRV  Data  (W 

(CS) 

Column  Name 

Format 

Description 

dt 

Timestamp 

Date/Time  of  measurement 

hrv_psn 

tinyint 

Pt  number  from  HRV  study 

hrv_evt 

integer 

Event  number  from  HRV  study 

q_admit_dt 

timestamp 

Admission  Date/Time 

T 

double 

Time  since  admission,  in  hours  &  decimal  fraction 

Table:  HRV.OXYl  -  HRV  Data  (W< 

CS) 

Column  Name 

Format 

Description 

dt 

Timestamp 

Date/Time  of  measurement 

hrv_psn 

tinyint 

Pt  number  from  HRV  study 

hrv_evt 

integer 

Event  number  from  HRV  study 

pao2 

integer 

Pa02 

saow 

integer 

SaC>2 

hct 

numeric(4,2) 

Hematocrit 

Table:  HRV.OXY2  - 

HRV  Data  (W< 

CS) 

Column  Name 

Format 

Description 

dt 

Timestamp 

Date/Time  of  measurement 

hrv_psn 

tinyint 

Pt  number  from  HRV  study 

hrv_evt 

integer 

Event  number  from  HRV  study 

svo2 

integer 

svo2 

o2_delivered 

integer 

02  Delivered 

o2_consumed 

integer 

02  Consumed 

dt 

Timestamp 

Date/Time  of  measurement 

hrv_psn 

tinyint 

Pt  number  from  HRV  study 

hrv_evt 

integer 

Event  number  from  HRV  study 

cobi 

integer 

Cardiac  Output  Bioimpedance 

cibi 

double 

Cardiac  Index  Bioimpedance 

hr_ekg 

numeric(4,l) 

HR  from  EKG 

map 

integer 

Mean  Arterial  Pressure 

ptc_02 

integer 

ptco2 

ptc_C02 

integer 

ptcco2 

FI02 

double 

FI02 

ptc02_FI02 

double 

ptco2/FI02 

sap02 

integer 

sapo2 

Table:  HRVJPT JLIST  -  HRV  Data  (WCS) 

Column  Name 

Format 

Description 

PSN 

numeric(7,0) 

Patient  Study  Number  (assigned  for  this  study) 

HRV.PSN 

tinyint 

Patient  number  from  HRV  study 

PF 

double 

Hospital  assigned  PF  number 

SURVIVED 

bit 

1  =  pt  survived,  0  =  pt  died 

AGE 

tinyint 

Age  in  years 

MALE 

bit 

1  =  male,  0  =  female 
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Column  Name 

Format 

Description 

HEIGHT 

Height 

WEIGHT 

SHI 

Weight 

DX 

919 

Diagnosis  &  notes 

Date/Time  of  admission 

Date/Time  of  operation  (if  any) 

DISCHARGE_DT 

timestamp 

BLOOD_LOS  S_EST 

smallint 

Blood  Loss  Estimate 

SSBEBBaB 

ARDS.DATE 

msmmm 

Date/Time  of  ARDS,  if  any 

BSA 

numeric(2,l) 

??? 

GCS  score 

Location 

Patient  location 

Notes 

Notes 

Table:  HRV_PT_LOCATION  -  HRV  Data  (WCS) 


Column  Name 

Format 

Description 

dt 

Date/Time 

hrv_psn 

Pt  number  from  HRV  study 

Location 

Pt  location 

Table:  icu_bp2  -  SICU  Data  (SICU) 


Column  Name 

Description 

emtek_id 

SICU  ID  number 

obj 

BP  source  (ART  or  BP) 

dia 

numeric(5,2) 

Diastolic  BP 

sys 

numeric(5,2) 

Systolic  BP 

mean 

numeric(5,2) 

Mean  Arterial  Pressure 

dt 

timestamp 

Date/Time  of  measurement 

Table:  icu_cardiac  -  S] 

1CU  Data  (SICU) 

Column  Name 

Format 

Description 

emtek_id 

integer 

SICU  ID  number 

svri 

numeric(5,l) 

Systemic  Vascular  Resistance,  Indexed 

svr 

numeric(6,2) 

Systemic  Vascular  Resistance 

CO 

numeric(5,2) 

Cardiac  Output 

ci 

numeric(5,2) 

Cardiac  Index 

dt 

timestamp 

Date/Time  of  measurement 

Table:  icu_demo  -  SICU  Data  (SICU) 


Column  Name 

Format 

Description 

PSN 

numeric(7,0) 

Patient  Study  Number  (assigned  for  this  study) 

HRVJPSN 

tinyint 

Patient  number  from  HRV  study 

PF 

numeric(8,0) 

Hospital  assigned  PF  number 

emtek_id 

integer 

SICU  ID  number 

admit_dt 

timestamp 

Admission  Date/Time 
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Column  Name 

Format 

Description 

icu_dt 

timestamp 

SICU  Admission  Date/Time 

ageadm 

smallint 

Age  at  admission 

bdate 

datetime 

Date  of  Birth 

weight 

numeric(5,2) 

Weight 

height 

numeric(5,2) 

Height 

gender 

char(l) 

Male  =  'M',  Female  =  'F 

survived 

char(l) 

Yes  =  ’Y',  No  =  ’N’ 

Table:  icujluids  -  SICU  Data  (SICl 

0 

Column  Name 

Format 

Description 

emtek_id 

integer 

SICU  ID  number 

vjn 

numerical) 

v_out 

numeric(8,l) 

Fluids  out  for  day 

v_net 

numeric(9,l) 

Net  fluids  for  day 

dt_day 

date 

Date  for  measurement 

Table:  icu_hr-  SICU  1 

Data  (SICU) 

Column  Name 

Format 

Description 

emtek_id 

SICU  ID  number 

hr 

Heart  Rate 

dt 

SSH 

Table:  icujabs  -  SICU  Data  (SICU) 

Column  Name 

Description 

emtek_id 

SICU  ID  number 

object 

Lab  test  abbreviation 

value 

flloat 

Lab  results 

dt 

timestamp 

Date/Time  of  measurement 

Table:  icu_rr  -  SICU  Data  (SICU) 


Column  Name 

■Kim 

Description 

emtekjd 

1  HI  1  "  III  — 

rr 

dt 

Table:  icutempr  -  SICU  Data  (SICU) 


Column  Name 

Format 

Description 

emtek_id 

tempr 

dt 

Table:toy_psn  (WCS,  SICU,  TR) 


Column  Name 

hezesbi 

Description 

psn 

HB 

pf 

numeric(7,0) 

2 


Column  Name 

Format 

Description 

hrv 

tinyint 

HRV  study  patient  number 

tr 

numeric(10,0) 

Trauma  Registry  pt  acct  number 

icu 

integer 

ICU  emtek_id  pt  number 

Table:  tr_demo-  TR  Data  (TR) 


Column  Name 

Format 

Description 

PSN 

numeric(7,0) 

Patient  Study  Number  for  this  study 

PF 

numeric(7,0) 

Hospital  assigned  PF  number 

HRVJPSN 

tinyint 

HRV  study  patient  number 

acctno 

decimal(20,5) 

Trauma  Registry  pt  acct  number 

emtek_id 

integer 

CU  emtek_id  pt  number 

er_enter_dt 

timestamp 

Date/Time  ER  entry 

er_exit_dt 

timestamp 

Date/Time  ER  exit 

er_los_calc 

numeric(6,2) 

er_exit_dt  -  er_enter_dt  in  hours 

er_los_rec 

numeric(5,2) 

ER  length  of  stay  as  recorded  in  TR 

rV_FLUIDS 

numeric(5,l) 

Total  IV  Fluids  from  ER 

Penetrating 

varchar(l) 

Penetrating  =  P,  Blunt  =  B 

surg_hrs_dt 

numeric(15,2) 

Hours  from  ER  admit  to  surgerv 

iss 

tinyint 

ISS 

ais_l_head 

tinyint 

AIS  Head  Score 

ais_2_head 

tinyint 

AIS  Face  Score 

ais_3_head 

tinyint 

AIS  Chest  Score 

ais_4_head 

tinyint 

AIS  Abdomen  Score 

ais_5_head 

tinyint 

AIS  Extremeties  Score 

ais_6_head 

tinyint 

AIS  External  Score 

sbp 

numeric(5,l) 

Systolic  BP 

dbp 

numeric(5,l) 

Diastolic  BP 

mbp 

numeric(5,l) 

Mean  BP 

Age 

numeric(5,2) 

Age 

Sex 

±ar(l) 

M/F 

Weight 

numeric(5,2)  ’ 

Weight  (kg) 

Survived  i 

:har(l) 

Y7N 

2.2.3  Main  Project  Database 

The  Main  Project  Database  contains  existing  medical  record  (EMR)  data  and  currently  consists  of 
689  trauma  patients  from  the  WCS  data  sets  and  corresponding  patient  data  from  the  TR.  As  of 
September  28,  2002,  the  SICU  database  has  been  sorted  and  pre-processed  for  inclusion  in  the  Main 
Database  but  had  not  yet  been  formally  included.  The  patients  studied  were  monitored  with  state-of-the- 
art  non-invasive  sensors  as  well  as  with  invasive  techniques  and  procedures.  Overall,  the  Database  is 
considered  to  be  very  high  in  quality  and  useful  in  prediction  because  the  diagnostic  measurements  are 
comprehensive,  and  the  prevalence  of  non-survivors  is  high.  The  Main  Project  Database  includes 
temporal  records  of  medications,  blood  component  therapies,  resuscitation  fluids,  and  procedures.  In 
addition,  the  location  of  the  patient  in  the  hospital  is  provided  as  a  function  of  time.  The  temporal  point 
of  reference  in  the  Database  is  the  time  of  admission. 
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Both  the  WCS  and  TR  database  involve  significant  numbers  of  keyed  entries  without  the  benefit 
of  a  supervising  program  to  prevent  unintended  characters/values  from  entering  the  database.  A  Visual 
Basic  program  was  written  to  examine  each  cell  of  the  WCS  Excel  sheet  for  correctness  of  data  format 
and  to  determine  whether  data  values  were  in  range.  This  procedure  helped  identify  errors  and  made  the 
transformation  of  the  Excel  Sheets  to  database  tables  easier.  Ultimately,  SQL  queries  of  the  data  sets 
within  the  Sybase  program  were  used  to  determine  whether  unreasonable  or  inconsistent  information 
was  contained  in  the  WCS  and  TR  data  sets.  In  the  case  of  the  WCS  data  set,  the  records  of  W.  C. 
Shoemaker  were  used  to  verify  the  data,  whereas  archived  patient  charts  were  examined  to  correct  the 
TR.  Correcting  the  TR  is  a  time  consuming  endeavor. 

In  Table  3,  we  characterize  the  data  from  the  WCS  and  TR  data  sets  that  currently  reside  in  the 
Main  Project  Database.  The  WCS  data  are  organized  into  several  individual  data  sets.  These  datasets 
have  been  used  in  several  focused  studies  conducted  in  the  past  at  LAC+USCMC  [e.g..  Refs.  6-11].  The 
percentage  of  patients  covered  by  various  physiological  parameters  is  presented  in  Table  3  to  help 
characterize  the  sensor  coverage.  Abbreviated  notations  for  the  physiological  variables  are  used  in  the 
tables  below;  detailed  definitions  of  the  parameters  are  provided  at  the  end  of  this  subsection  (pages  21- 
25). 


Table  4  summarizes  EMR  patient  data  processed  in  2002.  These  results  require  that  the  TR  be 
updated  by  hospital  personnel  before  the  new  data  can  be  fully  incorporated  into  the  Main  Project 
Database.  The  TR  updates  typically  lag  the  WCS  data  acquisition  by  -6-12  months.  During  the  current 
grant  year,  approximately  180  new  EMR  patient  histories  were  processed,  which  represents  an  increase 
of  -40%  in  the  number  of  trauma  patients  studied  compared  to  the  previous  year.  At  present,  the 
outcomes  of  some  of  the  patients  in  the  Surgical  ICU  5  data  set  are  not  known,  and  as  a  result  the  exact 
number  of  non-survivors  is  not  known. 

A  summary  of  physiological  parameters  versus  WCS  data  set  is  provided  in  Table  5  along  with  the 
percentage  of  patients  covered  by  each  measurement.  The  first  17  parameters  in  the  table  have  become 
standard  quantities  within  the  WCS  database  for  assessing  blood  flow  and  tissue  perfusion.  Several 
other  non-invasive  parameters  are  associated  with  the  heart  rate  variability  study.  The  standard  non- 
invasive  diagnostics  of  hemodynamics  include  electrical  bioimpedance  measurements  of  cardiac  output, 
pulse  oximetry,  and  transcutaneous  measurements  of  oxygen  and  carbon  dioxide  tension.  A  brief 
description  of  the  each  these  diagnostics  is  provided  below;  greater  detail  may  found  in  Refs.  6  and  12. 
The  Heart  Rate  Variability  Protocol  and  the  Surgical  ICU  2,  3,  and  4  studies  made  use  of  a  special 
purpose  instrument  that  measures  heart  rate  variability  parameters  (ANS-R1000,  Ansar  Inc., 
Philadelphia,  PA).  This  unit  yields  values  for  L(mean),  R(mean),  L/R(mean),  and  HR(mean)  shown  in 
Table  5. 

Bioimpedance  Measurements.  A  thoracic  bioelectric  impedance  device  (Yantagh,  Inc.,  Bristol,  PA) 
is  usually  applied  shortly  after  arrival  of  the  patient  in  the  emergency  department.  Noninvasive 
disposable  pre-wired  hydrogen  electrodes  are  positioned  on  the  skin,  and  three  EKG  leads  are  placed 
across  the  precordium  and  left  shoulder.  A  100  kHz,  4  mA  alternating  current  is  passed  through  the 
patient’s  thorax  by  the  outer  pairs  of  electrodes.  The  voltage  is  sensed  by  the  inner  pairs  of  electrodes, 
which  capture  the  baseline  impedance  (Z0),  the  first  derivative  of  the  impedance  waveform  (dZ/dt),  and 
the  EKG.  The  EKG  and  bioimpedance  signals  are  filtered  with  an  all-integer-coefficient  technology  to 
decrease  computation  and  signal  processing  time.  The  signal  processing  algorithm  uses  a  time- 
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Table  3.  W.  C.  Shoemaker  Data  Currently  in  the  Main  Database 


Study  Name:  Fluid  Resuscitation  Protocol 

Focus:  Effectiveness  of  resuscitative  fluids  on  tissue  perfusion. 

Number  of  Patients:  333,  78  non-survivors,  23.4%  non-survivors 
Earliest  Admission  Date:  7/27/96  (patient  1) 

Last  Admission  Date:  3/25/99  (patient  333) 

Recorded  Parameters  and  Percentage  of  Patients  Covered  by  Measurement 


m 

CIbi 

95.8 

HR 

95.8 

MAP 

94.3 

Sap02 

94.0 

Ptc02 

92.8 

PtcC02 

92.8 

FI02 

94.9 

Sa02 

59.2 

Sv02 

16.5 

HCT 

60.7 

D02I 

59.5 

V02I 

16.5 

Study  Name:  Heart  Rate  Variability  Protocol 

Focus:  Non-invasive  measurements  of  autonomic  nervous  system.  Heart  rate  variability  and 
respirator  activity. 

Number  of  Patients:  179, 43  non-survivors,  24.0%  non-survivors 
Earliest  Admission  Date:  8/1 8/99  (patient  2) 

Last  Admission  Date:  10/14/00  (patient  179) 

Recorded  Parameters  and  Percentage  of  Patients  Covered  by  Measurement 


COtd 

56.4 

COrb 

20.7 

CObi 

91.1 

CIbi 

91.1 

HR 

96.6 

MAP 

96.6 

Ptc02 

91.6 

PtcC02 

91.6 

FI02 

96.6 

tc02/F 

91.6 

Sap02 

96.6 

R(mean) 

98.9 

L/R(mean) 

99.4 

H 

1 

R(mean) 

00.0 

Pa02 

77.7 

Sa02 

77.7 

Sv02 

46.4 

HCT 

77.7 

ssjUl 

WSm 

V02I 

45.3 

Study  Name:  Surgical  ICU  2 

Focus:  Comparison  of  invasive  and  non-invasive  measurements  on  a  periodic  basis. 
Number  of  Patients:  67, 16  non-survivors,  23.9%  non-survivors 
Earliest  Admission  Date:  8/9/00  (patient  16) 

Last  Admission  Date:  5/31/01  (patient  67) 

Recorded  Parameters  and  Percentage  of  Patients  Covered  by  Measurement 


COtd 

67.2 

COrb 

25.4 

VC02 

25.4 

Pox 

25.4 

STAR  ' 
25.4 

rv 

55.2 

R1 

4C 

1.3 

CObi 

98.5 

CIbi 

98.5 

mmm 

tc02/F 

100.0 

Sap02 

98.5 

L(mean) 

11.9 

R(mean) 

11.9 

L/R(mean) 

11.9 

Pa02 

76.1 

IMEBMIili 

D02I 

55.2 

V02I 

52.2 

Qs/Qt 

52.2 

BE 

56.7 

PEEP 

50.7 

CVP 

38.8 

Study  Name:  Surgical  ICU  3 

Focus:  Comparison  of  invasive  and  non-invasive  measurements  on  a  periodic  basis. 
Number  of  Patients:  108,  21  non-survivors,  19.4%  non-survivors 
Earliest  Admission  Date:  3/11/01  (patient  1) 

Last  Admission  Date:  11/24/01  (patient  108) 

Recorded  Parameters  and  Percentage  of  Patients  Covered  by  Measurement 


CIbi 

HR 

MAP 

Sap02 

tc02/F 

Cild 

COrb 

W8M 

Pa02  82.4 

Sa02 

99.1 

100.0 

100.0 

100.0 

100.0 

46.3 

4.6 

82.4 

Sv02 

HCT 

D02I 

V02I 

Qs/Qt 

BE 

r 

I 

R(mean) 

L/R(mean) 

PtcC02 

49.1 

83.3 

50.0 

50.9 

50.0 

80.6 

17.6 

17.6 

17.6 

100.0 

Table  4.  W.  C.  Shoemaker  Data  Processed  in  Calendar  Year  2002 

(EMR  patient  data  processed  in  2002  require  that  the  TR  be  updated  before  the 
new  results  can  be  fully  incorporated  into  the  Main  Project  Database.) 

Study  Name:  Surgical  ICU  4 

Focus:  Comparison  of  invasive  and  non-invasive  measurements  on  a  periodic  basis. 

Number  of  Patients:  1 10,  28  non-survivors,  25.5%  non-survivors  (Preliminary) 

Earliest  Admission  Date:  11/26/01  (patient  1) 

Last  Admission  Date:  7/22/02  (patient  110) 

Recorded  Parameters  and  Percentage  of  Patients  Covered  by  Measurement 


CIbi 

100.0 

HR 

100.0 

m 

Sap02 

99.1 

PtcC02 

95.5 

tc02/F 

100.0 

Cltd 

47.3 

FI02 

100.0 

Pa02 

79.1 

Sa02 

79.1 

gji^| 

EHH 

HCT 

80.0 

D02I 

50.9 

V02I 

49.1 

Qs/Qt 

49.1 

BE 

75.5 

Study  Name:  Surgical  ICU  5  (in  Progress) 

Focus:  Comparison  of  invasive  and  non-invasive  measurements  on  a  periodic  basis. 

Number  of  Patients:  55,  at  least  9  non-survivors.  The  exact  number  of  non-survivors  cannot 
be  determined  as  of  10/09/02.  The  final  status  of  10  patients  is  not  known  at  present. 

Earliest  Admission  Date:  7/21/02  (patient  1) 

Lastest  Admission  Date:  9/28/02  (patient  55). 

Recorded  Parameters  and  Percentage  of  Patients  Covered  by  Measurement 


CIbi 

100.0 

HR 

100.0 

MAP 

100.0 

isa 

PtcC02 

100.0 

tc02/F 

100.0 

II 

FI02 

100.0 

Pa02 

73.3 

Sa02 

73.3 

Sv02 

40.0 

HCT 

73.3 

a 

V02I 

40.0 

Qs/Qt 

40.0 

BE 

70.0 

H 
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Table  5.  Summary  of  WCS  Physiological  Parameters  Versus  Data  Set  and 
the  Percentage  of  Patients  Covered  by  Each  Measurement 


Fluid  Heart  Rate 

Resus.  Variability  Surgical 

Parameter  Pr0t0C0|  protocol  ICU2 

_  07/96-03/99  08/99-10/00  08/00-05/01 


Cltd 

25.2 

Clbi 

95.8 

HR 

95.8 

MAP 

94.3 

EttJusSlHi 

94.0 

Ptc02 

92.8 

PtcC02 

92.8 

FI02 

94.9 

Sa02 

59.2 

Sv02 

16.5 

HCT 

60.7 

D02I 

59.5 

V02I 

16.5 

tc02/F 

92.8 

Pa02 


56.4 

67.2 

91.1 

98.5 

96.6 

100.0 

96.6 

100.0 

96.6 

98.5 

91.6 

92.5 

91.6 

92.5 

96.6 

100.0 

77.7 

76.1 

46.4 

52.2 

77.7 

76.6 

49.2 

55.2 

45.3 

52.2 

91.6 

100.0 

77.7 

76.1 

43.0 

52.2 

56.7 

20.7 

25.4 

20.7 

25.4 

91.1 

98.5 

97.8 

11.9 

98.9 

11.9 

99.4 

11.9 

100.0 

4.5 

VCQ2 

ETCQ2 

Pox 


Surgical 
ICU  3 

03/01  - 11/01 

Surgical 
ICU  4 

11/01  -07/02 

Surgical 

ICU  5 

07/02  -  09/02 
(In  Progress) 

46.3 

47.3 

36.7 

99.1 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

100.0 

99.1 

100.0 

100.0 

100.0 

100.0 

100.0 

95.5 

100.0 

100.0 

100.0 

100.0 

82.4 

79.1 

73.3 

49.1 

51.8 

40.0 

83.3 

80.0 

73.3 

50.0 

50.9 

40.0 

50.9 

49.1 

40.0 

100.0 

100.0 

100.0 

82.4 

79.1 

73.3 

50.0 

49.1 

40.0 

80.6 

75.5 

70.0 

frequency  distribution  (modified  Wigner  Distribution)  analysis  that  increases  the  signal-to-noise  ratio. 
Previous  studies  have  documented  satisfactory  correlation  between  thermodilution  and  bioimpedance 
cardiac  output  values  in  trauma  patients  [Ref.  13]. 

Pulse  Oximetry.  Routine  pulse  oximetry  (Nellcor,  Pleasanton,  CA)  is  used  to  continuously  assess 
arterial  oxygen  saturation  (Sap02).  Values  are  observed  and  recorded  at  the  exact  time  of  cardiac  index 
measurements.  Appreciable  or  sudden  changes  in  these  values  are  noted  and  confirmed  by  in  vitro 
arterial  oxygen  saturation  obtained  via  standard  blood  gas  analysis  [Ref.  13]. 

Transcutaneous  oxygen  and  carbon  dioxide  tensions.  Standard  transcutaneous  oxygen  tension 
(PtcP02)  measurements  (Novametrix  Medical  Systems,  Inc;  Wallingford,  CT)  are  made  continuously 
throughout  the  observation  period.  Values  are  noted  and  recorded  at  the  exact  times  of  cardiac  output 
measurements.  The  measurement  system  uses  the  same  Clark  polarographic  oxygen  electrode  routinely 
used  in  standard  blood  gas  analyses.  The  oxygen  tensions  are  determined  in  a  representative  area  of  the 
skin  surface  heated  to  44  C  to  increase  diffusion  of  oxygen  across  the  stratum  comeum  and  to  avoid 
vasoconstriction  in  the  local  area  of  the  skin  being  measured. 

Previous  studies  demonstrated  the  capacity  of  Ptc02  to  mirror  tissue  oxygen  tension  [Ref.  14]. 
Ptc02  has  been  shown  to  reflect  the  delivery  of  oxygen  to  the  local  area  of  skin;  it  also  parallels  the 
mixed  venous  oxygen  tension  except  under  late  or  terminal  conditions  where  peripheral  shunting  leads 
to  high  mixed  venous  hemoglobin  saturation  values  [Ref.  14].  Although  oxygen  tension  of  a  segment  of 
the  skin  does  not  replicate  the  state  of  oxygenation  of  all  tissues  and  organs,  it  has  the  advantage  of 
being  the  most  sensitive  early  warning  tissue  of  the  adrenomedullary  stress  response  and  the  adequacy 
of  oxygen  delivery  to  the  tissue.  Vasoconstriction  of  the  skin  is  an  early  stress  response  to  hypovolemia 
and  other  shock  syndromes. 

Transcutaneous  C02  (PtcC02)  monitoring  of  the  skin  surface  is  performed  with  the  standard 
Stowe-Severinghaus  electrode  ((Novametrix  Medical  Systems,  Inc;  Wallingford,  CT). 

Heart  rate  variability  measurements.  Heart  rate  (HR)  variability  and  respiratory  rate  (RR) 
variability  are  measured  by  spectral  analysis  techniques  with  the  ANS-R1000  instrument.  These 
parameters  are  indicative  of  autonomic  nervous  activity.  The  spectral  areas  of  variability  are  divided 
into  low  frequency  areas,  L(mean),  and  high  frequency  areas,  R(mean).  The  L(mean)  area  extends  from 
0.04  to  0.10  Hz.  This  area  reflects  primarily  the  tone  of  the  sympathetic  nervous  system  as  mediated  by 
the  cardiac  nerve.  R(mean),  also  referred  to  as  the  respiratory  area,  corresponds  to  a  0.12  Hz-wide 
frequency  interval  centered  on  the  fundamental  respiratory  frequency.  It  is  indicative  of  vagal  outflow 
reflecting  parasympathetic  nervous  system  activity.  The  L/R  ratio  reflects  the  balance  between  the 
sympathetic  and  parasympathetic  nervous  systems.  A  recent  study  [Ref.  8]  indicates  that  a  consistently 
positive  relationship  exists  between  HR  variability  and  L(mean)  during  sudden  surges  in  autonomic 
activity.  This  relationship  holds  to  a  lesser  degree  with  R(mean).  Heart  rate  variability  that  reflects 
autonomic  activity  is  associated  with  increased  MAP,  Cl,  and  HR,  but  decreased  tissue  perfusion  as 
indicated  by  the  Ptc02/FI02  ratio. 

The  decrease  in  the  heart  rate  variability  data  in  the  year  2002  (Surgical  ICU  4  and  5  data  sets)  is 
due  to  the  fact  that  the  ANS-R1000  instrument  is  currently  undergoing  hardware  and  firmware  upgrades 
at  Ansar  Inc.  The  future  role  of  this  instrument  in  WCS  studies  has  yet  to  be  determined. 
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Table  6  lists  medications,  blood  component  therapies,  resuscitation  fluids,  and  procedures  included 
in  the  WCS  database.  Patient  location  versus  time  can  also  be  determined  with  the  WCS  database. 
Possible  patient  locations  at  LAC+USCMC  are  presented  in  Table  7  below.  Finally,  the  TR  database 
parameters  selected  for  use  in  the  Main  Project  Database  are  listed  in  Table  8.  As  noted  earlier,  these 
parameters  provide  an  important  complement  to  the  WCS  and  SICU  databases. 


Table  6.  Listing  of  Medications,  Blood  Component  Therapies, 
Resuscitation  Fluids,  and  Procedures  in  the  W.  C.  Shoemaker  Data  Sets 


Medications 

Blood  Component  Therapies 
and  Resuscitation  Fluids 

Procedures 

DOP  (Dopamine) 

RBC  (Packed  Red  Blood  Cells) 

Intubation 

DOB  (Dobutamine) 

Plates  (Platelets) 

Entabation 

Morphine 

Cryo  (Cryoprecipitated  Anti¬ 
hemophilic  Factor) 

Dialysis 

NaHC03 

WBC  (White  Blood  Cells) 

Surgery  (start/closing/end 
times) 

Nitroprusside 

Colloids 

CPR 

T4  (Thyroxine) 

FFP  (Fresh  Frozen  Plasma) 

CT 

Nitroglycerin 

5%/25%  Albumin 

X-Ray 

Lasix 

Hespan 

Swan-Ganz  Catheter 

Insertion 

Atropine 

Crystalloids 

bagging  for  respirator 

25%  Mannitol 

Normal  Saline 

Chest  tube  insertion 

Pentobarbiturate 

Lactated  Ringer’s  Solution 

Angio  (Angiography) 

Suction 

Table  7.  The  Location  of  Patients  Versus  Time  is  Included  in  the 

W.  C.  Shoemaker  Data  Sets 


Locations 

_  Emergency  Department  (ED) _ 

_  Operating  Room  (OR) _ 

_ Radiology _ 

_ Surgical  Intensive  Care  Unit  (SICU) _ 

Continuously  Monitored  Area  (CMA)  — >  Step-Down  ICU 

Ward 
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Table  8.  Items  from  the  Trauma  Registry  that  are  in  the  Main  Database 


Injury  Date 

Abbreviated  Injury  Scale  (AIS) 

Injury  Time 

(l=Head,  2=Face,  3=Chest,  4=Abdomen,  Pelvis, 

Mechanism  of  Injury  Description 

5=Extremities,  6=Extemal  (e.g.,  skin)) 

Blunt/Penetrating 

Blood  (Total  blood/products  received  during  hospital  stay, 
including  Emergency  Department) 

Field  GCS 

Blood/Products 

Field  Airway  Type 

Autotransfuser 

Initial  Vital  Signs 

Total 

Time 

Discharge 

BP 

Discharge  Date 

RR 

Discharge  Time 

Assisted?  Y/N 

Prior  Phase 

HR 

Lived/Died 

Temp 

Organ  Donation  Y/N 

Weight 

Discharge  to: 

GCS 

Home 

Eye 

Other  Hospital 

Motor 

Trauma  Center 

Verbal 

Bum  Center 

Labs/Xray 

Rehab  Center 

C-Spine 

Skilled  Nursing  Facility 

CT  Head 

Morgue 

CXR 

Jail 

Pelvis 

AMA  (Patient  departed  Against  Medical 

Ultra  Sound 

Advice) 

CT  Abdomen 

Other 

CT  Spine 

Discharge  Diagnoses 

Facial  Series 

Complications  (Y/N) 

Other 

ARDS 

HCT 

Cardiac  Arrest 

ETOH 

Colonic  Anastomotic  Leak 

Toxicology  Serum 

Deep  Vein  Thrombosis  (DVT) 

Toxicology  Urine 

DVT  Lower  Extremity 

Procedure 

Disseminated  Intravascular  Coagulation  (DIC) 

ETT/CRIC/Trach 

Empyema 

ED  Thoracotomy 

Intra-abdominal  abscess 

DPL  (Peritoneal  Lavage) 

Jaundice  or  Hepatic  Failure 

Chest  Tube:  Lt/Rt 

Pancreatic  Fistula 

CPR  Duration 

Pneumonia 

IV  Fluids 

Pulmonary  Embolus  (PE) 

Pre-hospital 

Renal  Failure 

ED  IV  Fluids 

Septicemia 

Blood  Products 

Surgical  Dehiscence/Evisceration 

Autotrans 

Surgical  Wound  Infection 

None 

Other 
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Definitions  of  Physiological  Parameters  in  the 
W.  C.  Shoemaker  Data  Sets 


BE 

Base  Excess 
Blood  Gas  Lab  results 
(dimensionless,  pH) 

CIBI 

Cardiac  Index  Bio  Impedance 

IQ  machine 

l/min/mA2 

Integration  is  12  beats 

CITD 

Cardiac  Index  measured  by  Thermodilution 
Pulmonary  Artery  Catheter 
1  /min  /  mA2 

5  minutes  measurement  time 

COBI 

Cardiac  Output  Bio  Impedance 

IQ  machine 

1/min 

Integration  is  12  beats 

COrb 

Cardiac  Output  measured  by  C02  Rebreathing 
Non-Invasive  Cardiac  Output  (NICO) 

1/min 

COtd 

Cardiac  Output  measured  by  Thermodilution 

Swan-Ganz  TD  Catheter 

1/min 

~  5  minutes  integration  time 

CVP 

Central  Venous  Pressure 

Central  Venous  Line  -  average  pressure  from  central  line 
mm  Hg 

-  5  - 10  minutes  integration  time 
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D02I 

Oxygen  Delivery  Index 

Calculated  from  non-invasive  systems,  ICU  computer 

ml/minute/mA2 

15  minutes  measurement  time 

dZdt 

change  in  baseline 
IQ  machine 
Ohms/sec 

12-15  beat  average 

ETC02 

End  Tidal  C02 

NICO 

Torr 

~  1  minute  integration  time 

FI02 

Fractional  Inspired  Oxygen  Concentration 
Respirator 

Fraction  (1.00  =  100%  0.10  =  10%) 

Machine  setting,  not  a  measurement 

HCT 

Hematocrit 

Lab 

% 

HR 

Heart  Rate 

Bedside  EKG  if  no  dZdt  or  Z0  present,  otherwise  IQ  machine 
Beats  /  Minute 

10-15  cardiac  contractions  averaged 

HR_mean 

Heart  Rate 
HRV  machine 
beats/min 

15  minutes  integration  time 
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L_mean 

LFA  (Low  Frequency  Area) 

HR  spectrum  0.04  -  0. 10  Hz 
HRV  machine 
Hz 

15  minutes  integration  time 

LR_mean 

Ratio  of  L/R 
Unitless 

15  minutes  integration  time 

MAP 

Mean  Arterial  Pressure 

Bedside  monitor,  most  use  A-line,  some  use  cuff  (only  a  very  few). 
Not  indicated  in  data  set  which  was  used 
mm  Hg 

10-20  seconds  integration  time 

PA02 

Arterial  Blood  Gas  02  tension 

Lab  result 

Torr 

PAP 

Pulmonary  Arterial  Pressure 
PA  Swan  Ganz  Catheter 
mm  Hg 

-5-10  minutes  integration  time 

PEEP 

Positive  End  Expiratory  Pressure 

Respirator 

mm  Hg 

Pox 

(Sap02  from  NICO,  compared  with  Sap02  from  Colin) 

Pulse  Oximetry 
NICO 

% 

5-10  seconds  integration  time 
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PtcC02 

Transcutaneous  C02  Tension 

Novametrics 

Torr 

30  seconds  integration  time 

Ptc02 

Transcutaneous  Oxygen  Tension 

Novametrics 

Ton- 

30  seconds  integration  time 

Qs_Qt 

Physiological  Shunt  Pulmonary  Venous  Admixture 
amount  of  blood  going  through  lungs  w/o  being  oxygenated 
Calculated  from  Swan-Ganz 

% 

15  minutes  measurement  time 

R_mean 

HFA  (High  Frequency  Area) 

HR  spectral  curve  with  a  0. 10  Hz  window  around 
fundamental  respiratory  frequency 
HRV  machine 
Hz 

15  minutes  integration  time 

RR 

Respiration  Rate 

Respirator 

Breaths  /  Minute 

30  -  60  seconds  average 

SA02 

Saturation  of  Arterial  Hemoglobin 
Lab  result 

% 

SAP02 

Pulse  Ox 
Finger  Cuff 

% 

5-10  seconds  integration  time 
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STAR 

Goodness  of  ETC02  measurement 

NICO 

1-5  stars 

SV02 

Saturation  Venous  Oxygen 
Lab  Result 

% 

TC02/F 

Transcutaneous  02  /  FI02 
Calculated  from  TC02/FI02 
Torr 

30  seconds  integration  time 

TV 

Tidal  Volume 

Respirator 

liters 

30  -  60  seconds  integration  time 

V02I 

Oxygen  Consumption  Index 

Calculated  from  Swan-Ganz  TD  cath,  ICU  computer 

ml/minute/m2 

15  minutes  measurement  time 

zo 

Baseline  Bioimpedance 

IQ  machine 

Ohms 

12-15  beat  integration 
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2.2.4  New  Sensor  Data 

In  year  2  of  this  program,  EMR  measurements  made  with  the  CapnoProbe  sublingual  CO2 
monitoring  system  will  be  processed  and  added  to  the  database.  The  sublingual  approach  for 
measuring  CO2  was  pioneered  by  Max  Harry  Weil  (Institute  of  Critical  Medicine,  Palm  Springs, 
CA).  Subsequently,  Optical  Sensors  Inc.  (Minneapolis,  MN)  transformed  the  sublingual  probe  into 
an  FDA-approved  prototype  (Model  2000)  and  then  to  a  more  economical  handheld  instrument 
(M80).  This  latter  unit  is  marketed  and  sold  by  Nellcor  Puritan  Bennet,  Inc.  (Pleasanton,  CA). 

In  a  recent  clinical  study  involving  five  normal  human  volunteers  and  46  patients  with  acute 
life-threatening  illnesses  or  injuries,  sublingual  Pco2  Proved  t0  be  a  good  estimator  of  severity  of 
circulatory  shock  state  [15].  It  was  found  that  when  sublingual  PC02  exceeded  a  well-defined 
threshold  of  70  mm  Hg,  its  positive  predictive  value  for  the  presence  of  physical  signs  of  circulatory 
shock  was  1.00.  When  it  was  less  than  70  mm  Hg,  it  predicted  survival  with  a  predictive  value  of 
0.93.  A  majority  of  patients  in  whom  sublingual  Pco2  exceeded  70  mm  Hg  on  admission,  died  in 
the  hospital.  It  was  concluded  that  noninvasive  sublingual  Pcc>2  measurements  may  serve  as  a 
technically  simple  method  for  diagnosing  and  estimating  the  severity  of  circulatory  shock  states. 

As  part  of  the  current  project,  comparisons  will  be  made  between  transcutaneous  measurements 
of  CO2  and  sublingual  measurements  of  CO2.  The  potential  use  of  sublingual  CO2  as  an  outcome 
predictor  will  be  closely  examined. 

2.2.5  Improvements  In  Data  Exchange  at  LAC+USCMC 

During  the  current  grant  year,  the  automated  patient  monitoring  system  of  the  SICU  was 
enhanced  through  the  implementation  of  bi-directional  communications  with  the  hospital  pharmacy 
and  laboratory.  This  allows  the  complete  history  of  patient  medications,  fluids,  and  laboratory 
results  to  be  recalled  as  part  of  the  SICU  database.  Thus,  a  computerized  listing  is  now  available 
beginning  at  the  time  of  patient  admission,  and  this  information  can  readily  be  incorporated  into  the 
project  database. 

2.2.6  Database  Organization 

At  the  beginning  of  this  project  Geospace  Research,  Inc.  (GRI)  developed  a  custom  PC 
computer  architecture,  and  two  identical  PCs  were  designed  and  constructed.  One  was  placed  at 
GRI  and  the  other  was  located  at  LAC+USCMC.  These  computers  are  used  as  database  servers  for 
workstations.  Because  the  two  computers  are  identical,  the  mobile  database  disks  from  one 
machine  can  be  cloned  and  installed  on  the  other.  Thus,  updating  the  two  computer  systems  for  the 
latest  data  additions  is  a  simple  process. 

GRI  evaluated  two  database  programs  (Oracle  9i  and  Sybase  SQL  Anywhere  8.0)  before  finally 
selecting  the  Sybase  product  for  the  current  project.  There  are  many  reasons  why  Sybase  is 
preferable.  Sybase  has  a  much  better  tool  set  and  supports  a  publishable  database;  that  is,  it  allows 
portions  of  the  database  to  be  transferred  to  a  PDA,  a  laptop,  or  to  a  full  server.  Moreover,  Sybase 
requires  far  fewer  system  resources  than  Oracle  9i.  It  is  simple  to  administer  and  provides  all  of  the 
functions  necessary  for  the  project.  Sybase  8.0  is  a  workgroup  version  of  the  more  comprehensive 
Sybase  Enterprise  program.  The  advantage  of  Sybase  8.0  is  that  the  server  operates  peer-to-peer, 
and  a  server  operating  system  is  not  required.  If  necessary,  one  could  easily  upgrade  to  Sybase 
Enterprise  because  the  database  formats  are  identical.  However,  this  does  not  appear  necessary  on 
the  basis  of  our  first  year  results.  In  the  current  situation,  Windows  2000  is  used  as  the  database 
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server  operating  system  as  well  as  the  workstation  operating  system.  This  simplifies  computer 
operations.  Oracle  9i  (like  Sybase  Enterprise)  requires  a  server  operating  system.  For  our  particular 
application,  Oracle  9i  is  excessive  and  requires  a  steep  learning  curve.  For  comparison,  Oracle  9i  is 
distributed  on  14  CD’s,  whereas  Sybase  8.0  is  contained  on  a  single  CD.  Finally,  we  note  that  the 
SICU  database  discussed  above  is  administered  with  Sybase  Enterprise.  As  a  result,  the  SICU  data 
transfer  is  particularly  easy.  A  schematic  illustration  of  the  medical  database  currently  set  up  at 
GRI  is  shown  in  Figure  2  below. 

As  a  matter  of  policy,  ad  hoc  editing  of  data  in  the  database  is  not  permitted.  All  raw  data 
table  inputs  are  maintained  in  their  original  form.  When  necessary,  new  tables  containing 
transformed  and  edited  data  are  produced  in  addition  to  the  original  databases.  These  tables  are 
created  by  executing  SQL  scripts.  Desired  data  items  are  then  consolidated  into  a  single  database 
using  another  set  of  SQL  scripts.  This  policy  allows  our  working  database  to  be  recreated  at  any 
time,  with  all  editing  and  transformation  steps  clearly  documented  in  the  scripts.  Only  minor 
changes  in  the  SQL  scripts  are  necessary  to  run  them  under  another  database  engine. 

We  found  that  running  the  mining  and  analysis  tools  on  the  same  computer  was  problematic. 
(All  computers  discussed  here  have  a  Pentium  IV,  2-GHz  CPU  with  1-GB  memory;  the  operating 
system  is  Windows  2000  Professional.)  The  latest  MATLAB  software  (version  13.0)  now  has  a 
Java  user  interface.  The  MATLAB  Database  Toolbox  interfaces  with  Sybase  via  Microsoft's 
ODBC  drivers.  Running  Sybase,  Java,  MATLAB,  and  the  ODBC  on  the  same  machine  led  to 
unacceptable  processing  delays.  We  obtained  much  better  results  by  hosting  MATLAB  on  a  second 
machine  and  connecting  it  to  the  database  via  100  MBPS  Ethernet. 

At  present,  we  have  not  implemented  the  MATLAB  Web  Server  illustrated  in  Figure  2. 
Doing  so  in  the  future  would  allow  GRI  mining  library  routines  to  be  accessed  by  web  forms. 
Results  would  be  served  up  to  remote  users  either  as  plots  or  data  files.  Currently,  all  data  mining  is 
done  at  the  local  system  console  using  either  a  sequence  of  MATLAB  commands  or  MATLAB 
scripts  (i.e.,  m-files). 

2.2.7  Quality  Control 

In  principle,  the  construction  of  a  large  database  is  a  simple  matter,  but  in  practice  it  is  not. 
The  reason  for  this  is  that  errors  and  inconsistencies  in  the  incoming  data  must  be  isolated  and 
corrected.  As  the  database  becomes  very  large,  this  becomes  a  significant  problem.  The  two 
databases  that  are  most  susceptible  to  keying  errors  are  the  TR  and  WCS  database.  The  SICU  data 
is  less  problematic  because  the  data  acquisition  is  highly  automated,  a  review  is  required  for  data 
sign  off,  and  supervisory  software  is  used  to  prevent  erroneously  keyed  entries  from  being  inserted 
into  the  database. 

Because  the  WCS  database  is  in  the  form  of  Microsoft  Excel  spreadsheets,  a  series  of  data 
set  transformations  must  be  performed  to  convert  it  to  Sybase  format.  To  facilitate  this  process  and 
to  trap  errors,  each  cell  of  the  WCS  spreadsheet  is  initially  examined  with  the  aid  of  a  Visual  Basic 
program.  The  program  identifies  cells  that  are  improperly  formatted  (i.e.  as  a  text  string  instead  of  a 
numeric  value)  and  examines  all  data  to  make  sure  that  parameter  values  fall  within  reasonable 
limits.  If  a  data  value  is  suspect,  the  problem  is  addressed  using  the  original  data  records  of  W.  C. 
Shoemaker  to  resolve  the  issue.  Formatting  errors  prevent  the  translation  of  the  Excel  spreadsheet 
into  Sybase,  and  it  is  much  easier  to  identify  them  in  Excel  rather  than  through  terminating  errors  in 
the  translation  programs.  Excel  formatting  errors  are  often  subtle;  a  text  string  of  “10”  is 
indistinguishable  from  the  numerical  value  of  10  when  viewed  in  the  spreadsheet. 
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Figure  2.  Database  setup  for  the  LAC+USCMC  Medical  Trauma  Project. 
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The  TR  and  the  SICU  databases  are  supplied  as  relational  databases,  so  formatting  is  not  an 
issue.  However,  the  TR  database  structure  does  change  with  time.  Both  the  TR  database  and  the 
corrected  WCS  database  are  imported  into  the  database  program  in  the  form  of  tables.  Data  in  these 
“source  tables”  are  never  changed  or  edited.  Instead,  SQL  macros  are  used  in  conjunction  with  the 
tables  to  create  intermediate  data  tables.  When  new  TR  and/or  WCS  and/or  SICU  data  is  added  to 
the  system,  new  data  sets  are  generated,  and  analyses  are  performed  to  determine  whether  parameter 
values  are  reasonable  and  whether  overlapping  data  are  consistent.  For  example,  parameter 
distribution  functions  are  examined,  survivor/non-survivor  values  are  cross  referenced,  and  hard 
filtering  limits  are  placed  on  some  parameter  values  (e.g.,  blood  alcohol  level).  It  is  also  worth 
noting  that  a  patient  may  show  up  as  a  survivor  in  the  WCS  database,  and  a  non-survivor  in  the  TR. 
In  this  case,  the  patient  may  have  been  followed  for  two  weeks  for  the  WCS  database  and 
determined  to  be  survivor.  However,  on  occasion  patients  subsequently  undergo  sudden  downturns 
and  do  not  survive.  When  an  inconsistency  is  identified,  an  investigation  must  be  launched  to 
resolve  the  problem.  This  may  involve  locating  a  patient’s  archived  medical  chart,  which  can 
involve  a  search  lasting  several  days. 

Intermediate  databases  are  necessary  to  correct  errors  identified  in  the  original  databases  and 
to  normalize  certain  parameter  values.  For  example,  because  of  a  Y2K  problem  in  the  TR,  one 
patient’s  age  showed  up  as  137.  The  Y2K  algorithm  error  was  identified  and  corrected  as  part  of 
the  macro  that  created  an  intermediate  database.  The  original  data  remains  unchanged,  and  the 
macro  serves  as  documentation  of  the  changes  made.  The  macros  can  be  quite  long  because  they 
contain  information  about  specific  data  values  that  have  been  corrected,  data  transformations  that 
have  been  performed,  and  data  normalization  procedures  that  have  been  implemented.  In  the  latter 
case,  values  for  cardiac  output,  DO2,  VO2  etc.  may  have  to  be  indexed  prior  to  inclusion  in  the  main 
database.  Once  again  these  corrections  show  up  in  lines  of  SQL  code  in  a  macro. 

When  all  corrections  are  made  and  verified,  the  new  data  is  finally  incorporated  in  the  Main 
Project  Database.  Although  the  words  “seamlessly  add  new  patient  data  to  the  database”  were  used 
in  the  original  statement  of  work,  it  is  clear  that  the  necessary  quality  control  measures  prevent  this 
activity  from  being  truly  seamless  during  the  entry  process.  However,  the  final  Project  Database 
itself  is  seamless. 

2.2.8  Outcome  Prediction  With  the  Main  Project  Database 

Many  different  data  items  and  relationships  can  be  extracted  from  the  current  database. 
Figure  3  shows  a  plot  of  the  average  Cardiac  Index  (Cl),  Heart  Rate  (HR),  Mean  Arterial  Pressure 
(MAP),  percentage  oxyhemoglobin  saturation  (SapC>2),  and  transcutaneous  O2  tension  divided  by 
fractional  inspired  oxygen  concentration  (Ptc02/FI02)  versus  time  after  admission  to  the 
LAC+USCMC  emergency  department.  Patients  with  head  trauma  are  not  included  in  the  sample 
population  because  of  the  unique  type  of  hemodynamic  pattern  associated  with  head  injury. 
Standard  deviations  are  estimated  from  the  data  used  to  calculate  mean  values;  they  represent  the 
spread  in  the  distribution  of  data  values  about  the  mean.  Results  are  shown  in  Figure  3  for  survivors 
(blue)  and  non-survivors  (red).  A  good  predictor  yields  a  large  separation  between  the  parameter 
distribution  functions  of  survivors  and  non-survivors.  It  is  clear  that  PtcC>2/FI02  in  itself  is  a  fairly 
good  predictor  of  outcome,  particularly  at  the  earliest  times  after  the  trauma  victim  enters  the 
hospital.  The  large  Ptc02/FIC>2  standard  deviation  for  non-survivors  near  10  hours  is  caused  by  a 
low  patient  count  at  this  time.  Other  parameters  shown  in  Figure  3  are  not  good  predictors  because 
they  fail  to  produce  large  separations  in  survivor/non-survivor  distribution  functions.  As  noted  in 
Section  2.1.2  many  outcome  predictors  are  currently  being  assessed  as  part  of  a  more  generalized 
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Figure  3.  Mean  values  and  standard  deviations  of  five  hemodynamic  parameters  in  the 
Database.  Blue  designates  survivors,  whereas  red  represents  non-survivors.  Ptc02/FI02 
exhibits  the  greatest  population  separation  and  is  therefore  the  best  outcome  predictor. 
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study.  These  include  initial  heart  rate,  initial  MAP,  lowest  MAP,  initial  cardiac  index,  APACHE  II 
score,  predictions  of  discriminant  analysis,  probability  analysis  based  on  the  “nearest  neighbors”  of 
a  patient’s  “state  vector,”  PtcC^/FIC^,  and  the  combination  of  PtcC>2/FI02  and  PtcCC^.  At  present 
the  mathematical  “nearest  neighbor”  analysis  appears  to  be  the  best  technique  for  outcome 
prediction,  whereas  PtcC^/FIC^  has  emerged  as  the  best  single  parameter  predictor. 

Our  search  for  optimum  predictors  of  survival/non-survival  is  centered  on  finding 
physiological  parameters  that  yield  the  maximum  possible  separation  in  distribution  functions  of 
survivors  and  non-survivors.  In  this  regard,  one  examines  positive  (non-survival)  and  negative 
(survival)  predictability.  To  measure  the  effectiveness  of  the  predictor,  one  uses  indices  such  as 
sensitivity  (number  who  are  positive  and  test  positive/number  who  are  positive)  and  specificity 
(number  who  are  both  negative  and  test  negative/number  who  are  negative).  Initially  the  quantity  S 
=  [(1 -sensitivity)2  +  (1-  specificity)2]  is  minimized  to  establish  a  nominal  parameter  threshold  and 
two-sided  tests  are  used  to  determine  the  most  promising  predictors.  Subsequently,  trade-offs 
between  sensitivity  and  specificity  are  made  to  arrive  at  the  optimum  parameter  threshold.  This 
methodology  is  not  unlike  that  employed  to  interpret  laboratory  tests  or  to  diagnose  disease  with 
clinical  data.  In  our  case,  sensitivity  measures  how  well  the  physiological  quantity  identifies  those 
who  do  not  survive,  that  is,  how  sensitive  it  is.  If  a  test  has  a  high  sensitivity,  it  will  pick  up  nearly 
everyone  who  does  not  survive.  Specificity  measures  how  well  the  test  excludes  those  who  survive, 
that  is,  how  specific  it  is.  If  a  test  has  a  very  high  specificity,  it  won’t  misclassify  many  survivors  as 
non-survivors.  Typically,  one  analyzes  curves  of  sensitivity  versus  (1 -specificity)  to  establish  an 
optimum  cutoff  value. 

The  disadvantage  of  the  sensitivity  and  specificity  parameters  is  that  they  do  not  assess  the 
accuracy  of  the  test  in  a  clinically  useful  way.  However,  they  do  have  the  advantage  that  they  are 
not  affected  by  the  proportion  of  the  trauma  victims  that  do  not  survive,  which  is  called  prevalence. 
The  effect  of  lower  prevalence  is  much  as  one  would  expect:  the  lower  the  number  of  non¬ 
survivors  in  the  study,  the  more  certain  we  can  be  that  a  negative  test  indicates  survival  of  a  trauma 
patient,  and  the  less  sure  that  a  positive  result  really  indicates  that  the  trauma  victim  did  not  survive. 
In  the  current  investigation,  predictive  uncertainty  is  reflected  in  the  t  value  of  the  double-ended 
Student’s  t  distribution  and  in  the  fitting  errors  of  sensitivity  and  specificity  obtained  when  the 
quantity  S  is  minimized.  Because  the  prevalence  of  non-survivors  is  high  in  our  study  (20-25%), 
the  uncertainty  of  the  predictions  is  relatively  low. 

Software  has  been  written  to  automatically  examine  the  predictive  values  of  a  large  number 
of  parameters.  The  predictor  analysis  involves  least-squares  fitting  to  minimize  S  and  yields 
sensitivity  and  error  in  sensitivity,  specificity  and  error  in  specificity,  optimum  parameter  threshold, 
mean  of  parameter  value  in  non-survivors,  mean  of  parameter  value  in  survivors,  standard 
deviations  for  the  non-survivor  and  survivor  distribution  functions,  and  the  t  value  for  the  double 
ended  Student’s  t  Distribution.  The  main  W.  C.  Shoemaker  database  consisting  of  689  patients  is 
used  as  a  derivation  data  set  to  search  for  predictors.  A  separate  validation  data  set  is  being 
developed  to  confirm  the  predictors.  The  derivation  data  set  is  always  subject  to  bias  because  the 
same  data  used  to  derive  a  predictor  is  also  used  to  test  its  effectiveness.  Consequently,  there  is  a 
large  program  emphasis  on  securing  more  patient  case  histories.  This  is  evidenced  by  the  40% 
increase  in  total  number  of  EMR  trauma  patients  processed  during  the  current  grant  year. 
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3  Key  Research  Accomplishments 

1.  The  stage  1  development  of  a  large  national  trauma  database  for  studying  hemodyanamics 
including  blood  flow  and  tissue  perfusion  failure  was  completed. 

•  The  database  is  unbiased,  high  in  quality,  well  and  clearly  defined,  and  can  be 
independently  mined  by  members  of  the  combat-casualty-care  community. 

•  Mining  of  the  database  provides  insight  into  the  resuscitative  outcome  of  combat 
casualties  and  civilian  trauma  victims. 

•  The  database  can  be  used  to  help  determine  when  to  resuscitate,  how  to  resuscitate, 
and  what  the  end  point  of  resuscitation  should  be. 

2.  An  evaluation  of  outcome  predictors  indicates  that  a  mathematical  “nearest  neighbor” 
analysis  is  the  best  technique  identified  and  applied  to  the  database  for  predicting 
survival/non-survival,  whereas  Ptc02/FIC>2  has  emerged  as  the  best  single  parameter 
predictor. 

Outcome  predictors  greatly  facilitate  the  triage  process  for  combat  casualties  and 
measure  the  effectiveness  of  the  resuscitation  effort. 

Outcome  predictors  together  with  the  analyses  of  hemodynamic  patterns  allow 
hospitals  to  identify  patients  who  are  at  risk  in  the  earliest  stage  of  the  therapeutic 
process,  adjust  therapies  to  improve  outcomes,  and  promptly  determine  whether  the 
new  therapy  will  lead  to  survival. 

4  Reportable  Outcomes 
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1.  Shoemaker  WC,  Wo  CCJ,  Botnan  A,  Bayard  DS,  Jelliffe  RW:  Development  of  a 
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outcome.  IEEE  Transactions  on  Automatic  Control  2001;  36:  No.  9, 231-236 

2.  Bayard  DS,  Botnan  A,  Shoemaker  WC,  Jelliffe  RW:  Stochastic  analysis  of  therapeutic 
modalities  using  a  database  of  patient  responses.  IEEE  Transactions  on  Automatic  Control 
2001;  36:No.  9, 439-444 
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patients  by  noninvasive  hemodynamic  monitoring.  Chest  2001;  120:528-537 

4.  Olinski  M,  Shoemaker  WC,  Reis  E,  Kerstein  M:  Current  controversies  in  shock  and 
resuscitation.  Surg  Clin  North  America  2001;  81:1217-12^2 

4 

5.  Shoemaker  WC:  New  approaches  to  trauma  management  using  severity  of  illness  and 
outcome  prediction  based  on  noninvasive  hemodynamic  monitoring.  Surg  Clin  N  Am  2002; 
82:245-255 

6.  Kern  J,  Shoemaker  WC:  Meta-analysis  of  hemodynamic  optimization  in  high  risk  patients. 
Crit  Care  Med  2002;  30: 1686-1692 

7.  Shoemaker  WC,  Bayard  DS,  Wo  CCJ,  Botnan  A,  Chan  L,  Jelliffe  RW,  Djuth  FT,  Belzberg 
H:  Noninvasive  hemodynamic  monitoring  of  emergency  patients  for  outcome  prediction  by 
a  stochastic  control  program.  Crit  Care  Med  2002,  Submitted 

8.  Fathizedeh  P,  Shoemaker  WC,  Wo  CCJ,  Colombo  J:  Autonomic  activity  in  trauma  patients 
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4.2  Presentations 


1.  Shoemaker  WC:  Noninvasive  monitoring  of  high  risk  patients  using  a  new  stochastic 
analysis  and  control  system.  Aukland,  New  Zealand  April  8, 2002. 

2.  Shoemaker  WC:  Noninvasive  monitoring  of  high  risk  patients  using  a  new  stochastic 
analysis  and  control  system.  Christchurch,  New  Zealand.  April  10, 2002. 

3.  Shoemaker  WC:  Noninvasive  monitoring  of  high  risk  patients  using  a  new  stochastic 
analysis  and  control  system.  Wellington,  New  Zealand,  April  12, 2002. 

4.  Shoemaker  WC:  Noninvasive  monitoring  of  high  risk  trauma  patients  using  a  new  outcome 
predictor,  stochastic  analysis  and  therapeutic  decision  support  programs.  Sacramento  County 
Trauma  Service,  Roseland  Hospital,  Sacramento  CA.  May  30, 2002. 

5.  Shoemaker  WC:  Noninvasive  monitoring  of  high  risk  patients  using  a  new  stochastic 
analysis  and  control  system  for  outcome  prediction  and  therapeutic  decision  support  system. 
Frankfort,  Germany,  June  6,  2002. 

6.  Shoemaker  WC:  Outcome  prediction  and  decision  support  program  in  high  risk  trauma  and 
surgical  patients  using  a  new  stochastic  analysis  with  a  control  system  and  a  therapeutic 
decision  support  program.  Berlin  University,  Germany,  June  10, 2002. 

7.  Shoemaker  WC:  Noninvasive  hemodynamic  monitoring  of  high  risk  patients  using  a  new 
stochastic  analysis  for  outcome  prediction  therapeutic  decisions.  Mexican  Pediatric  Society, 
International  Conference,  Mexico  City,  Mexico,  June  20, 2002. 

8.  Djuth  FT,  Belzberg  H,  Shoemaker  WC,  Elder  JH,  Zhu  J,  Wo  CCJ:  An  open  trauma  database 
for  studying  hemodynamics  including  blood  Flow  and  tissue  perfusion  failure,  ATACCC 
Meeting,  St.  Pete  Beach,  FL,  September  9-13, 2002. 

9.  Belzberg  H,  Elder  JH,  Djuth  FT,  Oder  D,  Shoemaker  WC:  The  developing  of  a  Trauma 
Outcome  Data  Analysis  (TODA)  tool:  answering  old  questions  with  new  techniques: 
ATACCC  Meeting,  St.  Pete  Beach,  FL,  September  9-13, 2002. 

10.  Shoemaker  WC,  Belzberg  H,  Bayard  DS,  Wo  CCJ,  Botnen  A,  Djuth  FT,  Choi  A,  Jelliffe, 
RW:  Noninvasive  hemodynamic  monitoring  of  trauma  patients  for  outcome  prediction  and 
decision  support  by  a  stochastic  control  analysis  and  display,  ATACCC  Meeting,  St.  Pete 
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4.3  Databases 

1.  EMR  test  (Toy)  database  consisting  of  53  trauma  patients  common  to  the  WCS  database 
(Heart  Rate  Variability  data  set),  the  Surgical  Intensive  Care  Unit  database,  and  the  Trauma 
Registry  database. 

2.  EMR  Main  Project  Database  consisting  of  689  trauma  patients  from  the  WCS  data  sets  and 
corresponding  patient  data  from  the  Trauma  Registry.  The  patients  were  monitored  with 
state-of-the-art  non-invasive  sensors  as  well  as  with  invasive  techniques  and  procedures. 
This  database  is  considered  to  be  very  high  in  quality  because  the  diagnostic  measurements 
are  comprehensive,  and  the  prevalence  of  non-survivors  is  high. 
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5  Conclusions 

The  evolving  trauma  database  developed  as  part  of  this  project  serves  as  a  tool  for 
understanding  the  resuscitation  process.  The  further  development  of  appropriate  methodologies  and 
tools  for  mining  the  database  will  greatly  improve  the  quality  of  trauma  care  for  both  the  combat 
casualty  and  the  civilian  trauma  victim.  The  organized  database  supports  analyses,  evaluations,  and 
data  mining  by  diagnostic  and  covariate  categories,  as  well  as  by  temporal  patterns  of  hemodynamic 
patterns  in  survivors  and  nonsurvivors.  Our  preliminary  results  indicate  that  the  new  outcome 
predictor(s)  are  greater  than  90%  correct  within  the  first  12  to  24  hours  after  emergency  department 
admission.  On  the  basis  of  such  measurements,  a  therapeutic  decision  support  and  predictive 
modeling  system  is  being  developed  that  relies  primarily  on  the  hemodynamic  patterns  of  survivors 
and  their  response  to  therapy. 

5.1  Suggested  Changes 

During  the  first  year  of  work,  it  became  clear  that  more  existing  medical  record  (EMR) 
patient  data  will  greatly  enhance  the  trauma  database  and  support  complex  methodologies  used  to 
mine  the  database.  In  response  to  this  need,  additional  EMR  data  were  processed  during  the  first 
grant  year.  The  concept  of  prospective  data  acquisition  suggested  in  the  Year  2  statement  of  work 
will  not  generate  the  large  number  of  trauma  cases  necessary  to  meet  the  emerging  needs  of  the 
program.  However,  the  processing  of  an  ever  increasing  supply  of  EMR  patient  data  will. 
Therefore,  we  suggest  that  paragraph  two  of  the  Year  2  statement  of  work  be  replaced  by  the 
following: 

“During  the  second  grant  year,  additional,  newly  acquired,  existing  medical  record  patient 
data  will  be  processed,  analyzed,  interpreted,  and  incorporated  into  the  Main  Project  Database 
without  the  possibility  of  patient  identification.” 

A  related  e-mail  communication  from  Dr.  Maryann  F.  Pranulis,  AMEX/AMEDD,  is 
included  as  an  attachment  to  this  annual  report. 

5.2  Importance 

The  development  of  the  current  trauma  database  is  expected  to  greatly  improve  the  quality 
of  trauma  care  for  both  the  combat  casualty  and  the  civilian  trauma  victim.  Taken  as  a  whole,  this 
new  database  and  associated  mining  tools  allow  optimum  resuscitation  approaches  to  be  developed 
and  provide  a  means  for  measuring,  validating  and  implementing  new  physiological  sensor 
strategies.  At  present,  it  also  appears  that  a  new  therapeutic  decision  support  system  for  trauma 
physicians  will  be  an  important  spin-off  of  this  project.  The  Main  Project  Database  is  open  to  other 
investigators,  contains  an  ever  increasing  number  of  high-quality  well-documented  trauma  cases, 
and  accommodates  new  state-of-the-art,  FDA-approved  sensors  and  techniques.  As  a  result,  the 
time  line  required  to  develop  more  effective  therapies  and  translate  this  research  into  clinical 
practice  will  be  much  shorter  than  in  previous  studies  of  this  nature. 
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7  Appendices 

1.  Recent  e-mail  communication  from  Dr.  Maryann  F.  Pranulis,  human  subjects  protection 
scientist,  AMEX/AMEDD,  concerning  the  approval  of  exempt  status  and  the  elimination  of  a 
prospective  study. 

2.  Copies  of  eight  publications  listed  in  Section  4  above. 
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Date:  Wed,  02  Oct  2002  15:53:09  -0400 

From:  "Pranulis,  Maryann  F  Dr  AMDEX"  <Maryann.Pranulis@DET.AMEDD.ARMY.MIL> 
Subject:  RE:  A-10847,  DAMD1 7-01 -2-0823,  proposal  00089002, 

Studies  of  Tissue  Perfusion  Failure  at  LAC  &  USCM... 

To:  'Howard  Belzberg'  <belzberg@usc.edu> 

Cc:  "Stotler,  Karen  S  Ms  USAMRAA"  <Karen.Stotler@DET.AMEDD.ARMY.MIL>, 

"Vandre,  Robert  H  COL  USAMRMC"  <Robert.Vandre@DET.AMEDD.ARMY.MIL>, 
"Zadinsky,  Julie  K  COL  USAMRMC"  <Julie.Zadinsky@DET.AMEDD.ARMY.MIL>, 
"Pranulis,  Maryann  F  Dr  AMDEX"  <Maryann. Pranulis@DET.AMEDD.ARMY.MIL> 

X-Mailer:  Internet  Mail  Service  (5.5.2656.59) 

Importance:  high 

Dr.  Belzberg: 

It  was  a  pleasure  to  talk  with  you  today  and  to  learn  that  the  retrospective  component  of 
this  award  has  yielded  sufficient  data  to  delay  or  outright  drop  doing  a  prospective  study. 

As  I  understand  from  our  discussion,  you  have  made  no  changes  to  the  retrospective 
component  that  was  approved  as  exempt  on  29  September  2001 .  So  long  as  there  are  no 
changes  in  venue,  investigators,  procedures,  or  data  base  that  will  be  used,  and  you  are 
still  using  de-identified,  existing  medical  record  data,  the  previous  approval  of  exempt 
status  is  still  valid.  However,  it  should  be  noted  that  the  cited  approval  restricted  you  from 
recruiting  and/or  enrolling  any  human  subjects  or  using  any  prospectively  collected  human 
subjects  data.  The  approval  was  also  only  for  use  of  the  trauma  data  base  at  the  Los 
Angeles  County  and  University  of  Southern  California  Medical  Center  Surgical  Intensive 
Care  Unit.  These  restrictions  are  still  in  effect. 

If  there  are  any  modifications  in  the  previously  approved-as-exempt  protocol,  then  you  are 
required  to  either  submit,  for  HSRRB  review  and  approval,  a  new  protocol  or  a  request  for 
an  amendment  to  the  approved  version  prior  to  implementing  the  change.  If  the  USC-  IRB 
changes  the  designation  from  exempt  to  non-exempt,  please  notify  this  office  and  you  will 
be  advised  of  the  procedures  that  need  to  be  followed. 

Best  wishes  for  continued  success  in  this  important  endeavor.  If  you  have  any  questions  or 
concerns  about  the  human  subjects  aspects  of  this  project,  please  let  me  know. 

Maryann  F.  Pranulis,  RN,  DNSc 
Human  Subjects  Protection  Scientist 
AMDEX 

maryann.pranulis@det.amedd.army.mil 

(301)619-6240 

FAX  (301)  619-7803 

Message-ID:  <5.0.2.1. 1.2002100210471 1.00ab7140@email.usc.edu> 

From:  Howard  Belzberg  <belzberg@usc.edu> 

To:  maryann.pranulis@det.amedd.army.mil 
Subject: 

Date:  Wed,  2  Oct  2002  14:06:58  -0400 
MIME-Version:  1 .0 

X-Mailer:  Internet  Mail  Service  (5.5.2656.59) 


Content-Type:  multipart/alternative; 

boundary=" — _=_NextPart_003_01  C26A4D.551 3C420" 

Dear  Dr.  Pranulis, 

In  regard  to  our  project  (contract  no.  DAMD  17-01-2-0070)  we  plan  to  modify  our  plans  for 
year  2.  After  review  of  our  first  year's  acquisition  of  patients  using  retrospective  de- 
identified  data,  and  discussions  with  our  sponsor,  we  do  not  plan  to  persue  the  prospective 
component  of  the  project  at  the  present  time.  These  changes  will  be  outlined  in  our  annual 
report  to  be  submitted  by  October  28, 2002.  In  light  of  these  changes  we  would  request 
that  the  previously  granted  exemption  be  continued. 

Thank  you, 

Howard  Belzberg,  M.D. 


Development  of  a  Hemodynamic  Database  in  Severe  Trauma  Patients 
to  Define  Optimal  Goals  and  Predict  Outcome. 

WC  Shoemaker  ",  CCJ  Wo0,  A  Botncn2’,  DS  Bayard",  RW  Jclliffc2, 

Department  of  Surgery  and  Anesthesia  ",  Laboratory  of  Applied  Pharmocokinetics  ", 
University  of  Southern  California.  School  of  Medicine.  Los  Angeles  CA  90033  and 
Jet  Propulsion  Laboratory  ",  Pasadena  CA  91109 


Abstract 

Noninvasive  hemodynamic  monitoring  systems  provide  continuously  monitored  on  line 
displays  of  data  from  emergency  department  admission  to  the  OR.  and  to  the  1CIJ  for  early 
recognition  of  circulatory  dysfunction  in  acute  emergency  conditions.  The  net  cumulative 
defats  of  cardiac  index  are  estimated  by  thoracic  electric  bioimpedance,  arterial  hypoxemia 
ts  measured  by  pulse  oximetry,  and  tissue  perfusion  is  reflected  by  transcutaneous  p()r  Based 
on  a  large  database,  survival  was  satisfactorily  predicted  by  discriminant  analysis  and  bv  a 
new  stochastic  analysis  and  control  program. 


1.  Introduction 

We  studied  sequential  hemodynamic  patterns  of  surviving  and  nonsurviving  postoperative 
patients,  and  used  the  temporal  circulatory  patterns  of  survivors  in  the  first  P  hours 
postopcrativcly  alter  severe  trauma  as  a  first  approximation  to  define  optimal  goals  ofdicrapv 

and  inc  pauems  o!  nonsurvi  vors  as  early  warning  signs  t>i  potcniial  disasters  1 1 

Goal-directed  therapy  has  been  used  as  a  strategy  lor  hieh-risk  surucrv  and  acute  illness 
such  as  trauma  |l-7|.  Allliough  most  shock  studies  have  demonstrated  that  earlv  diagnosis  and 
vigorous  therapeutic  interventions  were  keys  to  success,  early  diagnosis '  of  drculatorv 
dyslunetion  is  difficult;  because  traditionally  shwk  has  been  rceoenized  bv  imprecise  siern 
and  observer-dependant  subjective  symptoms.  Shock  is  easy  to  rccouni/e  in  its  later  siauesbv 
hypotension,  oliguria,  and  collapse,  but  by  that  time,  therapv  is  often 'incITeciivc 
_  Noninvasi  ve  hemodynamic  monitoring  technology  has 'now  become  available  throughout 
the  hospital  including  the  emergency  department  (ED),  operating  room  (OR),  posi-uncsihcsia 
recovery  room.  ICU,  step  down  units,  hospital  floors,  prehospital  areas,  and  doctors*  offices 
Noninvasive  monitoring  is  easier,  faster,  safer,  cheaper  and  equally  sensitive.  18-14 1. 

,  htajor  lesson  learned  Iroin  earlier  studies  is  that  effective  resuscitation  of  acute  life- 
threatening  emergencies  requires  achieving  optimal  physiological  goals  as  early  as  possible:  if 
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circulatory  mechanisms  can  be  identified  earlier  and  treated  more  vigorously  to  specified 

physiologic  criteria,  outcomes  should  be  improved. 

We  have  developed  a  large  database  of  invasive  and  noninvasivc  hemodynamic  »h>»«  onng 
of  acutely  ill  trauma  and  surgical  emergencies  beginning  in  the  emergency  department  (ED):  a) 
to  describe  early  survivor  and  nonsurvivor  patterns  of  emergency  patients  in  icnns  or  cardiac, 
pulmonary,  and  tissue  perfusion  deficiencies;  b)  to  measure  quantitatively  die  net  cumu  alive 
amount  S'  deficit  or  excess  of  the  monitored  functions  that  correlate  with  survival  or  death 
and  c)  to  use  discriminant  analysis  to  predict  outcome  and  evaluate  the  brolog.cal  significance 

of  monitored  deficits  ( 1 3). 

2.  Methods 

2.1  Development  of  databases 

Databases  for  acute  injury  patients  (N=316).  elective  high  risk  surgery,  noninvasivcly 
monitored  patients  (N=374)  and  invasivcly  monitored  patients  (N=809),  septic  shock  patients 
(N=378),  and  hemorrhage  only  patients  (N=  87)  have  been  develop  to  describe  the ^pnniary 
injuries  or  illnesses,  the  hemodynamic  patterns  by  invasive  and  nonnwas.vc  nicll« od .  J md  tl^cn 
Hitcomcs,  including  survival/nonsUrvival,  organ  lailurc,  other  complications,  tune  ol  ED 
admission  hospital  days,  1CU  days  and  estimated  hospital  costs.  Noninvasivc  monitoring  was 
begun  in  the  ED  and  followed  to  die  OR.  radiology  department  and  1CU;  invasive  pufiiKinaiy 
artery  catheterization  (PAG)  was  instituted  when  clinically  indicated  after  arrival  in  the  ICU, 

have  developed  a  laree  database  wid.  M  10.000  do*  lines  descrihtne  die  clnneal  and 
hemodynamic  patterns  of  30  clinical  subsets  with  specified  co-morbid  conditions  18 
hemodynamic  variables,  and  over  2500  therapeutic  interventions  with  values  obtained  bclorc. 
during  and  after  each  therapy  given  onc-al-a-limc,  except  for  clinical  exigencies. 

The  sequential  patterns  were  described  in  time  elapsed  from  onset  or  illness  or  '» 

ED  admSion.  Weighting  criteria  to  specify  the  nearest  neighbor,  were 

diagnosis  co-morbid  conditions,  and  hcmodynam.c  patterns  ol  survivors  vs  nonsurv.vors.  I  he 
p  valucs  ^  diffcrenccs  in  the  temporal  patterns  of  survivors'  and  nonsurvivors  values  were 
Lsed  as  weighting  criteria  in  each  subset,  for  example,  for  subsets  of  patients  with  and  without 
head  injuries,  blunt  vs.  penetrating  trauma,  truncal  and  nontruncal  trauma.  age  stratifications, 
nrior  cardiac,  respiratory,  hepatic,  and  renal  dysfunction  or  organ  lailurc.  etc.  ... 

P  Finally  the  hemodynamic  responses  to  standardized  test  doses  of  whole  blood  or  packed  red 
cell  SSs  ^en  over  a  specified  period  of  lime  (usually  .-hour,  provide  quantitative 
measures  or  each  patient’s  cardiac  reserve  capacity.  This  is  expressed  as  changes  in  cardiac 
index  relative  to  corresponding  change  in  PA  occlusion  (wedge,  pressure  or  central  venous 
pressure  (Starling’s  myocardial  performance  curve)  before  and  alter  a  blood  translusion  given 
over  a  one-hour  period.  Criteria  for  satisfactory  vs.  limited  cardiac  functional  reserve  capacity 
have  been  developed  lor  each  clinical  subset  and  each  time  period. 

2.2  Noninvasive  Hemodynamic  Monitoring  Systems 

The  monitoring  systems  arc  an  improved  noninvasivc  thoracic  electric  bioimpedanee 
cardiac  output  device,  a  pulse  oximeter  to  measure  arterial  Itcmoglob.n  saturation,  a 
noninvasive^blood  pressure  device,  and  transcutaneous  Oj  and  C02  sensors  to  relied  tissue 
ncrfusion/oxygenalion.  These  continuously  monitored  noninvasivc  measurements  were  used  to 
circulatory  patterns  in  151  consecutively  monitored  severely  injured 
pmients  beginning  with  admission  to  the  emergency  department  in  a  university  run  county 
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hospital.  The  net  cumulative  deficit  or  excess  of  each  monitored  parameter  was  calculated  as 
die  cumulative  difference  from  die  normal  value  versus  the  time-integrated  monitored  curve 
for  each  patient.  The  deficits  Of  cardiac,  pulmonary,  and  tissue  perfusion  functions  were 
analyzed  in  relation  to  outcome  by  discriminant  analysis  and  crossvalidated.  * 

Noninvasivc  monitoring  compared  favorably  with  invasive  thcrmodilulion  catheter  f 
monitoring  in  high  risk  patients  with  surgery,  blunt  trauma,  gunshot  wounds,  head  injuries, 
sepsis,  strokes,  drug  overdose,  myocardial  infarction,  and  acute  gastrointestinal  bleeding. 
Differences  between  impedance  and  thcrmodilulion  cardiac  output  estimations  were  more  than 
offset  by  the  continuous  on-line  display  of  data  that  allowed  instant  recognition  of  changes  in 
the  course  of  a  patient’s  illness,  calculation  or  die  deficit  or  each  monitored  variable,  and 
evaluation  of  therapeutic  responses. 

3.  Results 

3.1  Therapeutic  Goals 

The  specific  monitored  goals  of  therapy  based  on  empirical  findings  of  survivors  were: 
cardiac  index  >  4.5  L/min/m2;  systolic  blood  pressure  >  120  niniHg;  pulse  oximetry  >  96%. 
transcutaneous  oxygen  PtcCVFiOi  ratio  >200;  heart  rale  <  100  bcats/inin  (1,3). 

3.2  Survivor  and  Nonsurvivor  Hemodynamic  Patterns  in  High-risk  Surgery,  1  raunin, 
and  Sepsis 

In  nonsurvivors,  the  initial  hemodynamic  changes  consisted  of  reduced  cardiac  index, 
oxygen  delivery,  and  oxygen  consumption.  Shortly  after  surgery,  the  survivors  had 
hemodynamic  and  oxygen  transport  values  greater  than  normal,  and  greater  than  those  in  the 
nonsurvivors.  Similarly,  in  severely  traumatized  patients  and  in  both  medical  and  surgical 
septic  shock  patients,  survivors  also  had  increased  hemodynamic  and  oxygen  transport  values, 
suggesting  (hat  these  responses  were  compensating  lor  prior  inadequacies  of  tissue  perfusion 
and  oxygenation. 

3.3  Quantitative  Assessment  of  Continuously  Monitored  Noninvasive  Variables  as  Net 
Cumulative  Amount  of  Excess  or  Deficit 

Many  recordings  show  considerable  variability  which  may  obscure  the  underlying  pattern. 

To  overcome  this  problem  we  evaluated  the  net  cumulative  deficit  or  excess  in  cardiac1 
pulmonary,  and  tissue  perfusion  function  data  by  integrating  over  lime  the  area  between  the 
monitored  curve  and  the  normal  or  optimal  values.  In  nonsurviving  patients,  the  mean  values 
as  well  as  the  net  cumulative  deficits  were  significantly  greater  for  cardiac  index,  tissue 
perfusion,  and  arterial  saturation  than  in  survivors,  even  tliough  the  survivors'  anil 
nonsurvivors'  MAP  values  were  not  significantly  different. 

3.4  Outcome  Prediction  by  Discriminant  Analysis 

The  data  were  analyzed  by  discriminant  function  (15).  The  following  variables  were  selected 
by  (lie  stepwise  discriminant  analysis  (PROC  STEPDISC):  a)  cumulative  PtcO^/FiOj,  b) 
Glasgow  coma  score,  c)  cumulative  Sap02  values  and  d)  cumulative  cardiac  index.  Based  on 
the  classification  function  generated  for  each  of  these  four  variables  in  PROC  DISCRIM, 
discriminant  function,  Z,  was  derived  for  an  individual  patient:  Z  =  0.001  la  +  0.33001’ 
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.  0  0656c  +  0  0423d;  where  “a”  represents  cumulative  PteCtyFiOj  value,  “b”  represents 

represents  Lndadve  SapO, Z 
cardiac  index  value.  Classification  of  the  survivors  was:  Z>2.36.  Half  of  the :  eases  were 
randomly  selected  to  serve  as  the  calibration  data  set;  the  other  half  was  used  as  the  validation 
Z2ZS2ZZZ- '•  above,  there  were  l«  (18.4*)  ,„iselass,r, canons  ,,,  (he 
cross-validation  study,  compared  with  23/151  (15.2%)  of  the  series  as  a  whole.  This  was 
considered  to  be  in  reasonably  satisfactory  agreement  (14).  Using  stepwise  discnininant 
analysis;  95  %  of  the  survivors  and  62%  of  the  nonsurvivors  were  corrccUy  classified  in  the 
early  period  after  the  initial  resuscitation. 


Figure  1.  Classification  Summary  for  the  Series  (N=151) 

4.  Conclusion 

Hemodvnamic  data  obtained  by  pulmonary  artery  (PA)  thcrmodiluiion  (Swan-Ganz) 
catheter  have  shown  significantly  higher  hemodynamic  and  oxygen  transport  values  in  survivors 
compared  with  nonsurvivors  (1,2.4-12).  Classically,  shock  is 

livDOtcnsion,  and  collapse  in  the  late  stage,  when  therapy  is  least  effective.  However  initial  y 
Jrfrccoenizcd  shock  by  subjectively  evaluated  signs  and  symptoms  such  as  skin  color,  skin 

weak  Uirrady  pulse.  and  vial  sign,  Tire  convenuonal  lher*euue 

!£S£s  ton  toconccl  (hcan  superlielal  nranifaaadons  . -  8«=y  an:  reeugnwed.  bn. 

« rrcrn» 

suDPlv  information  previously  available  only  in  die  cadnacrizalion  lab  under  sienlc  opcralmg 
SSSS  o, by  L  PA (Swan-oana,  oati^er -Jr ^““tTusS 

be  used  Ibrouglrout  Ore  bosptial  in  dre 
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early  period  immediately  after  emergency  department  admission.  They  provide  continuous  on¬ 
line  displays  rather  than ’’snapshots”  at  infrequent  intervals. 

Most  importantly,  in  the  period  during  and  immediately  after  resuscitation,  noninvasivc 
hemodynamic  patterns  or  patients  with  severe  trauma  arc  closely  related  to  outcome.  We  have 
used  these  patterns  to  develop  outcome  predictors  by  discriminate  analysis  to  elucidate 
underlying  mechanisms  (13).  Also  we  now  propose  to  supplement  or  replace  discriminate 
analysis  With  a  stochastic  control  program,  which  is  a  more  powerful  research  tool  applicable 
to  iiKirc  extensive  use  of  muniiMcd  variables  in  a  wider  variety  of  acutely  ill  patients. 

The  data  show  that  the  initial  hemodynamic  findings  are  low  now  or  unevenly  distributed 
microcirculatory  How  witli  poor  tissue  perfusion/oxygenation;  these  are  initially  precipitated 
by  many  factors  including  hypoxemia,  hypovolemia,  acidosis,  and  the  adrcnomcdullary  stress 
response.  The  poorly  perfused,  acidolic  capillary  endothelial  wall  activates  macrophages, 
stimulates  the  systemic  immune  response  syndrome  (SIRS),  and  produces  reactive  oxygen 
substances  (ROS),  that  may  contribute  to  die  development  of  ARDS  and  other  organ  failures.  _ 
Our  hypothesis  is  that  early  identification  of  the  initiating  hemodynamic  mechanisms  allows 
their  correction  with  concomitant  reversal  of  shock  and  improved  outcome.  Preliminary  data 
demonstrated  that  the  stochastic  control  program  provides  objective,  quantitative  information 
on  the  likelihood  of  survival  or  death  before  and  alter  each  therapy.  In  essence,  the  stochastic 
control  program  provides  early  warning  or  initiating  mechanisms  at  times  when  they  can  be 
reversed,  evaluates  the  relative  efficacy  of  specific  therapy  to  reverse  early  deleterious 
uicchanisms.  and  thereby  provides  the  basis  lor  understanding  of  the  contributory  role  ol  these 

hcnuKlynamic  mechanisms.  ■ 

This  analysis  is  based  on  the  assumption  that  it  may  be  easier  and  more  efficacious  to 
identify  initiating  hemodynamic  mechanisms  of  shock  and  lethal  organ  failure  by  evaluation  ol 
early  hemodynamic  patterns  and  their  transitions  alter  therapeutic  interventions.  Transitional 
states  alter  therapeutic  interventions  also  may  be  evaluated  in  terms  or  possible  mechanisms; 
for  example,  cardiac  reserve  capacity  may  be  inferred  from  changes  in  myocardial 
performance  after  a  standardized  transfusion  of  I  unit  of  whole  blood  given  over  a  l-hour 
period  or  a  standardized  test  dose  of  dobulaniinc.  A  distinction  is  made  between  initialing 
hemodynamic  mechanisms  and  later  immunochemical  or  other  biochemical  mediators  of  organ 

failure.  .  ,  , 

Outcome  predictors  will  be  used  to  fine  tune  evolving  hemodynamic  patterns  and  to  observe 

the  effects  or  modifying  potentially 'lethal  mechanisms  by  therapy  early  in  the  course  of  acute 
illness.  The  assumption  is  that  there  arc  different  types  and  degrees  ol  hemodynamic 
dysfunction  that  need  objective  measures  to  track  their  presence  and  intensity.  An  objective 
outcome  measure  representing  the  probability  ol  survival  should  be  a  useful  means  to  evaluate 
the  evolution  of  acute  illnesses  that  lead  to  organ  failures.  Trauma  and  high-risk  surgical 
patients  were  selected  because  of  their  varying  hemodynamic  patterns,  severity  of  illness,  and 
associated  co-morbid  conditions. 

5.  References 

1.  Shoemaker  WO,  Appel  PI.  Krnm  IIB.  ct  at:  Prospective  trial  of  supranormal  values  of  survivors  as 
therapeutic  goals  in  high-risk  surgical  palicnls.  Clicsl  1988;  94:1 176-86 

2.  Boyd  O.  Grounds  M,  Bcnncll  I):  Preoperative  increase  ol  oxygen  delivery  reduces  mortality  in  high  risk 
surgical  patients.  JAMA,  1993;  270:  2699 

3.  Shoemaker  WC.  Appel  PI-  Kram  IIB:  Hemodynamic  and  oxygen  transport  responses  in  survivors  and 

nonsurvivprs  ol  high  risk  surgery.  Crit  Care  Med  1993;21:977-90 


BEST  AVAILABLE  COPY 


236 


4.  Bishop  Mil,  Shoemaker  WC,  Kram  MB,  et  aJ;  Prospective  trial  of  survivors  values  of  cardiac  index,  oxygen 
delivery,  and  oxygen  consumption  as  resuscitation  endpoints  in  severe  trauma.  J  Trauma  1995;  38:780-87 

5.  Wilson  J,  Woods  I,  Fawcett  J,  et  al:  Reducing  the  risk  of  major  surgery:  randomized  controlled  trial  of 
preoptimization  of  oxygen  delivery.  Brit  Med  J  1999;318:1099-1 103 

6.  Lobo  SMA,  Salgado  PF,  Castillo  VGT,  cl  al:  HfTccts  of  maximizing  oxygen  delivery  on  morbidity  and 
mortality  in  high  risk  surgical  patients.  Crit  Care  Med  2000;  28:  3396-3404 

7  Boyd  O:  The  high  risk  surgical  patient  -  where  are  we  now?  Clin  Intensive  Care  2000;  Special  issue;  3- 10 

8.  Wang  X,  Sun  H,  Adamson  D,  et  al:  An  impedance  cardiography  system:  A  new  design.  Ann  Biomed  Hug 
1989;  17:535-56 

9.  Shoemaker  WC,  BeJzberg  ,  Wo  CCJ,  Miltzman  D,  et  al:  Multiccntcr  trial  of  noninvasivc  monitoring  systems 
as  alternatives  to  invasive  monitoring  of  emergency  patients.  Chest  1998;!  14:  643-52 

10.  Raaijmakcrs  E,  Faes  TJC,  Scholton  RJPM,  et  al:  A  mcla-analysis  of  three  decades  of  validating  thoracic 
impedance  cardiography.  Crit  Care  Med  1999;  27:1203-13 

1 1.  Tremper  KK,  Waxman  K,  Shoemaker  WC:  Effects  of  hypoxia  and  shock  on  transcutaneous  POj  values  in 
dogs.  Crit  Care  Med  1979;  7: 526 

12.  _  Tremper  KK,  Iluxtable  RF:  Dermal  heat  transport  analysis  for  transcutaneous  02  measurements.  Acta  Ancslli 

Scand  1978  (suppl);  68: 4 

13.  Shoemaker  WC,  Wo  CCJ,  Chan  L,  Ramiconc  15,  et  al:  Outcome  prediction  of  emergency  patients  by 
noninvasive  hemodynamic  monitoring.  Chest  In  Press 

14.  Bayard  DS,  Jelliffe  RW,  Botnen  A,  Wo  CCJ,  Shoemaker  WC:  Stochastic  Analysis  and  Control  using  a 
database  of  patient  responses  to  therapy. 


BEST- AVAILABLE  COPY 


Stochastic  Analysis  of  Therapeutic  Modalities  Using  a  Database  of 

Patient  Responses 

D  S  Bayard,  A  Botncn,  W  C  Shoemaker,  R  Jclliffc 

Laboratory  of  Applied  Pharmacokinetics.  USC  School  of  Medicine.  Los  Angeles  CA 

I 

Abstract 

„  °  ne"'  me,hoJf',r  anal ysis  and  control  which  Joes  no,  reunite 

/  II -I  and  methods  of  dynamic  programming  for  stochastic  control  given  in  /.*  //■/  /.  '  ^ 

1.  Concept  of  State 

si^liliEissts 

given  patient  (c «  cardinr  imw  wt  i  c  tvpcs  of  measurements  obtained  on  a 

tensions  etc  »  ThrnT  k  hkxxl  Prcssurc-  Pulse  oximetry.  transcutaneous  O.  and  CO 

ddlnC  *hc  vector  af a 

first  integral  j  y,a, .  as  follows  SCC,’ml  l"”c  dcrivalivcs  .VJ,  •  and  with  i,s 


•^)  =  |  y.^Xy,(r*  ly'M  *|  vr/f......y,(f,  iy,  jjy, J  (1) 

2.  Implicit  Dynamic  model 

.  »>K"k  TO«i,»  „f  ,hc  pal™-,  ,,  „  linlc  ,  .  his 

m”li”ar  ,1>"™ica,  srs*cm  *«*  »;l 

■vt*l  a  f\Xt.Ut,p.lVt)  ,,, 
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For  simplicity,  p  is  assumed  to  be  drawn  from  a  finite  set  formed  by  enumerating  all  useful 
combinations  of  co  van  ales, 

pe  {p,,...,  /?„}  (3) 

Both  the  parameter  vector  p  of  covariatcs  and  process  noise  u-  help  to  explain  the  variability  of 
responses  seen  in  the  database.  Specifically,  the  covariatcs  help  to  distinguish  gross  differences  in 
responses  due  to  major  categories  of  complications  (c.g.,  with  vs',  without  head  injuries,  blunt  vs. 
penetrating  trauma,  truncal  vs.  nomruucul  trauma,  age  stratifications,  etc.)  On  the  other  hand, 
process  noise  helps  to  explain  smaller  differences  between  patients  with  identical  covariatcs  but 
with  different  responses  to  tire  same  therapy. 

3.  Control  Inputs 

It  will  be  assumed  that  there  is  only  a  finite  set  of  M  control  inputs  in  (2)  that  can  be  applied 
to  the  system.  Specifically,  the  control  u\tk  j  at  timet*  is  assumed  to  be  drawn  from  the  finite  set, 

u(tk)e  ,  (4) 

where  i/f  arc  controls  (therapeutic  nmdafitics)  such  as  fluid  therapy  Willi  crystalloids,  with 
colloids,  with  whole  bkxxl  or  packed  RBC’s,  or  various  vasoactive  drugs  such  as  dobutaminc, 
dopamine,  and  other  possible  modalities. 

4.  Nearest  Neighbors 

Once  the  state  vector  and  covariatcs  arc  defined,  me  key  concept  of  a  "nearest  neighbor”  can 
be  put  forth.  Given  the  state  a*  and  a  certain  co variate  p ,  the  N  nearest-neighbor  states 
(denoted  as  «  n (A* ,pj)  arc  defined  as  the  states  X  that  arc  closest  to  xk  in  the  database 

and  which  share  the  same  covariate  vector  p .  Here,  a  measure  of  “closeness”  is  conveniently 
defined  in  terms  of  the  quadratic  distance  between  die  state  and  its  neighbors, 

cl (.v.  A* )  =  (a*  -  xY  W (.v*  -  a)  (5) 

where  vv  is  an  appropriately  chosen  weighting  matrix. 

5.  Performance  Measure 

The  performance  measure  is  the  probability  of  survival.  For  a  patient  in  a  given  state  x  with 
covariate  vector  p ,  the  survival  probability  is  denoted  by  S(x,p).  It  is  evaluated  simply,  by 
extracting  the  N  nearest  neighbor  states  to  x  of  patients  with  the  same  covariate  vector  p  and 
by  noting  the  fraction  of  them  that  survived.  For  example,  if  l\}  of  the  N  nearest-neighbor 
states  to  x  survived,  then  the  survival  probability  is  given  as, 

S(x,p)  =  -ij-  (6) 
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6.  Applications 


6.1  Data  mining 

One  can  now  trace  through  individual  patient  histories  in  the  database  and  compute  the 
associated  probability  of  survival  i>[x,p)AS  a  function  of  time.  A  large  change,  either  up  or 
down,  indicates  a  potentially  significant  role  of  the  associated  therapeutic  modality,  and  its 
underlying  physiologic  mechanisms.  'I  he  entire  database  can  be  mined  lor  therapeutic  modalities 
having  the  most  eff  icacy  or  harm  under  specific  conditions  using  this  technique. 

6.2.  Real-Time  Diagnostic  Tool 

The  stochastic  analysis  is  not  restricted  to  subjects  contained  within  the  database.  II  can  also  bi¬ 
ased  on  any  real  patient  during  his/her  clinical  care.  In  this  case,  the  survival  probability  S{x,u) 
is-  computed  by  tracking  (he  subject's  nearest  neighbors  in  the  database  and  providing  to  the 
clinician  a  real-time  feedback  signal  which  indicates  the  condition  of  the  patient  as  a  function  ol 
lime,  and  the  efficacy  of  any  therapeutic  modalities  administered,  or  to  be  considered  at  that  time. 

6.3. Stochastic  Control  Optimal  Therapy 

By  classifying  a  subject's  nearest  neighbors  by  the  therapeutic  modalities  employed,  the 
survival  probabilities  , p )  can  be  computed  separately  for  each  modality  us.  A  useful 

method  for  stochastic  control  then  follows  by  administering  the  therapy  associated  with  the 
highest  value  or  F{u,txktp).  Intuitively,  this  therapeutic  modality  is  optimal  (and  a  worthwhile 
choice  in  practice)  in  tlic  sense  dial  it  has  the  highest  survival  probability  among  the  patient  s 
nearest  neighbors  in  the  database. 

The  stochastic  analysis  and  control  approach  is  summarized  in  Figure  I.  Examples  ol  the 
second  application  above  (i.c.,  use  as  a  real-time  diagnostic  tool)  will  be  given  in  the  next  section. 


«*(*»•  f) 

Figure!.  Slochaslic  analysis  and  control  synthesized  using  a  database  of  responses. 
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7.  Examples 

In  this  section,  two  examples  arc  given  demonstrating  the  use  of  the  stochastic  analysis  tool  for 
rtal-time  diagnosis.  For  this  purpose,  a  hemodynamic  database  or  responses  is  developed  from 
<hc  records  of  316  acutely  ill  patients  15 1. 

The  covariate  associated  with  truncal  injury  is  chosen.  The  stale  is  chosen  as  a  vector  whose 

dements  ate  composed  of  the  time  t  elapsed  since  admission,  and  observations  y,,l  -  1, . 7 

and  their  first  and  second  derivatives-  v, ,  y, ,  and  first  integral  j  y,ai  for  each  or  the  following 
observations:  Cl  (cardiac  index);  HR  (heart  rale).  MAP  (mean  arterial  pressure),  pulse  oximetry 
arterial  oxygen  saturation  Sap02 ),  transcutaneous  oxygen  Ptc07 ;  transcutacneous  carbon 
dioxide  PtcC07 :  and  HOT  (hematocrit).  The  resulting  slate  has  29  dimensions. 

Data  sets  of  two  patients  were  removed  from  the  database  and  treated  as  new  patients.  The 
probability  of  survival  (bottom  traces)  is  calculated  using  the  stochastic  analysis  approach  and 
plotted  as  a  function  of  time  in  Figure  2  for  patient  1  and  in  Figure  3  for  patient  2. 


Sgure2.  Patient  1  Data  -  Upper  row:  Cardiac  output  (Cl);  Second  Row:  Mean  arterial  pressure 
IAP);  Third  row:  Pulse  oximetry  (Sap02);  Fourth  row:  Transcutaneous  FI02  ratio;  Lowest  row: 
Survival  Probability  computed  tor  this  patient  by  the  stochastic  analysis  and  control  program. 
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Figure  3  Patient  2  Data  -  Upper  row;  Cardiac  output  (Cl);  Second  Row:  Mean  arterial  pressure 
(MAP);  Third  row:  Pulse  oximetry  (Sap02);  Fourth  row:  Transcutaneous  FI02  ratio;  Lowest  row: 
Survival  Probability  computed  for  this  patient  by  the  stochastic  analysis  and  control  program} 

tn  cadi  figure,  the  dark  ureas  arc  regions  of  deficit,  where  the  patients*  measurements ‘Tell 
below  set  thresholds  into  non-ideal  territory.  The  more  time  spent  in  such  regions  (i.c.,  building 
up  an  “integrated  deficit”),  the  more  reason  lor  concern  ior  the  well-being  o(  die  patient.  This 
integrated  deficit  statistic  has  been  studied  separately  and  has  emerged  as  a  strong  indicator  of 
patient  survivability  in  earlier  studies  |6).  It  is  used  here  to  compare  with  the  new  survival 
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probability  statistic.  In  cadi  figure,  tinic  in  hours  from  ED  admission  arc  noted  on  the  dark 
horizontal  line  below  die  survival  probability  row.  The  start  and  end  of  the  surgical  operation 
(OP)  are  noted.  Each  therapy  is  noted  in  the  oblong  boxes,  and  times  of  the  beginning  and  end  of 
each  therapy  arc  noted  by  die  vertical  dotted  lines. 

In  the  data  of  patient  1,  (Figure  2)  continuing  deficits  in  Cl  and  the  transcutaneous  OVFIO2 
ratio  resulted  in  a  progressive  lowering  of  this  patient’s  survival  probability  beginning  5  hours 
after  admission.  The  patient  subsequently  died  on  the  9'”  postop  day.  The  survival  probability 
correctly  indicates  this  degrading  trend,  becoming  smaller  with  lime  and  ending  with  a  survival 
probability  of  less  than  50%. 

In  the  data  of  patient  2  (Figure  3)  there  arc  relatively  small  deficits  in  Cl  values,  which  were 
usually  above  die  normal  value  of  3  L/inin/m2,  but  not  at  the  optimal  values  of  4  L/min/m2.  Tlicrc 
was  also  slight  hypotension,  despite  a  relatively  normal  Cl  and  the  relatively  low  tissue  perfusion 
indicated  by  the  PtcOVFIO2.  Despite  these  mild  to  moderate  deficiencies,  this  patient  survived. 
The  survival  probability  correctly  predicts  this  optimistic  trend.  Following  an  initial  low  it 
remained  above  75%  toward  the  middle  and  at  the  end. 

In  both  examples,  it  is  seen  that  the  survival  probability  has  provided  a  reliable  real-time 
diagnostic  signal  to  indicate  the  condition  of  die  patient. 
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Objectives:  We  used  noninvasive  hemodynamic  monitoring  in  the  Initial  resuscitation  beginning  in  the 
emergency  deportment  (ED)  for  the  following  reasons:  (1)  to  describe  early  survivor  and  nonsurrivor 
patterns  of  emergency  patients  hi  terms  of  cardiac,  pulmonary,  and  tissue  perfusion  deficiencies;  (2)  to 
measure  quantitatively  the  net  cumulative  amount  of  deficit  or  excess  of  the  monitored  functions  that 
correlate  with  survival  or  death;  and  (3)  to  explore  the  use  of  discriminant  analysis  to  predict  outcome  and 
evaluate  tiie  biological  significance  of  monitored  deficits. 

Methods:  This  is  a  descriptive  study  of  the  feasibility  of  noninvasive  monitoring  of  patients  with  acute 
emergency  conditions  in  the  ED  to  evaluate  and  quantify  hemodynamic  deficits  as  early  as  possible.  The 
noninvasive  monitoring  systems  consisted  of  a  bloJmpe dance  method  for  estimating  cardiac  output 
together  with  pulse  oximetry  to  reflect  pulmonary  function,  transcutaneous  oxygen  tension  to  reflect  tissue 
perfarion,  and  BP  to  reflect  the  overall  circulatory  status-  These  continuously  momtored  noninvasive 
measurements  were  used  to  prospectively  evaluate  Circulatory  patterns  in  151  consecutively  monitored 
severely  injured  patients  beginning  with  admission  to  the  ED  in  a  umversity^run  county  hospital.  The  net 
cumulative  deficit  or  excess  of  each  monitored  parameter  was  calculated  as  the  cumulative  difference  from 
the  normal  value  vs  the  time-integrated  monitored  curve  for  each  patient  The  deficits  of  cardiac, 
pulmonary,  and  tissue  perfusion  functions  were  analyzed  in  relation  to  outcome  by  discriminant  analysis 
and  were  cross-validated. 

RgbuUh:  The  mean  (X  SEM)  net  cumulative  excesses  (+)  or  deficits  (-)  from  normal  in  surviving  vs 
nonsurviving  patients,  respectively,  were  ns  follows:  cardiac  index  (Cl),  +81  ±  52  vs  —232  ±  J.38  L/m 
(p  =  0.037);  arterial  hemoglobin  saturation,  -X  ±  0.3  vs  -8  ±  2.6%/h  (p  =  0.006):  and  tissue  perfusion, 
+313  ±  83  vs  -793  2:  175,  mm  Hg/h  (p  =  0.001 ).  The  cumulative  mean  arterial  BP  deficit  for  survivors  was 
—  10  2=  13  mm  Hg/h7  and  for  nonsurvivors  it  was  —57  ±  24  Dim  Hg/h  (p  -  0.078). 

Conclusion*:  Noninvasive  monitoring  systems  provided  continuously  monitored  on-line  displays  of  data  in 
the  early  postadmission  period  from  the  ED  to  tlx?  operating  room  and  to  the  ICIT  for  early  recognition  of 
circulatory  dysfunction  m  short-term  emergency  conditions.  Survival  was  predicted  by  discriminant 
analysis  models  based  on  the  quantitative  assessment  of  the  net  cumulative  deficits  of  Cl»  arterial 
hypoxemia,  and  tissue  perfusion*  which  were  significantly  greater  in  the  nonsurvivors. 

(GUEST  2001;  12Q:52?i-537) 

Key  words:  hemodynamic  monitoring:  multicomponent  noninvasive  circulatory  monitoring:  outcome  prediction’;  p.»lse 
oximetry,  tamport]  hemodynamic  patterns;  transcutaneous  oygan  tension 

Abbreviations:  Cl  =  cardiac  index:  ED  =  emergency  department;  Fl0a  “  fraction  of  inspired  oxygen;  GCS  =  Glasgow 
com 9  scale;  MAP  =  mean  arterial  BP;  OR  *=  operating  room;  PAC  *■=  pulmnuaiy  artery  catheter;  Sa02  **  arterial  oxygen 
saturation;  icPoa  tmnscutfincouS  oxygon  tension 
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Hemodynamic  bedside  monitoring  by  pulmonary 
artery  catheters  (FACs)  has  been  considered  by 
many  as  the  "gold  standard*’  for  critically  ill  patients, 
but  its  usefulness  has  been  challenged, 1,  7  particu¬ 
larly  .in  the  late  stages  of  illness  after  the  onset  of 
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organ  failures.  A  meta-analysis  by  Boyd  and  Hayes* 
showed  no  outcome  improvements  in  seven  random¬ 
ized  studies  of  patients  who  entered  the  ICU  after 
organ  failure  or  sepsis  had  occurred,  but  there  was 
significantly  improved  outcome  in  six  other  random¬ 
ized  studies,  plus  two  more  recent  studies,9'10  when 
PAC-directed  therapy  was  given  early  or  prophylac- 
tically.  Since  time  may  be  important  in  the  initial 
resuscitation  and  management  of  emergency  pa¬ 
tients,  noninvasive  monitoring  is  proposed  as  an 
alternative  approach  to  identify  and  correct  hemo¬ 
dynamic  deficiencies  at  the  earliest  possible  time. 
Previous  studies  have  documented  satisfactory  cor¬ 
relation  between  thermodilution  and  bioimpedance 
cardiac  output  values  for  trauma  patients  in  the 
emergency  department  (ED),  the  operating  room 
(OB),  and  the  ICU,  The  mean  (±  SEM)  bias  and 
precision  in  the  ED  were  -0.055  X  0.75  I/min/m2.11 

In  the  present  study,  we  monitored  severely  in¬ 
jured  emergency  patients,  beginning  in  the  ED  and 
continuing  in  the  radiology  department,  the  OR,  and 
men  in  the  ICU.  Acute  injury  was  studied  because 
time  factors  are  important  and  the  time  course  of 
circulatory  events  could  be  monitored  ftom  the  time 
of  hospital  admission.1112  Continuous  visual  displays 
of  monitored  data  were  used  to  evaluate  rapidly 
changing  patterns  during  unstable  emergency  con¬ 
ditions.  Second,  we  time-integrated  the  differences 
between  the  monitored  curve  and  normal  values  or 
reference  values  reflecting  “optimal"  goals  derived 
from  the  patterns  observed  throughout  the  time 
course  of  previous  series  of  survivors  of  acute  severe 
illnesses  or  operations.13-21  We  then  calculated  the 
net  cumulative  excesses  or  deficits  of  each  monitored 
variable  for  each  patient  and  for  the  survivors  and 
nonsurvivors.  Finally,  we  explored  the  use  of  dis¬ 
criminant  analysis  to  predict  outcome  based,  on  these 
calculated  cumulative  deficits. 


Materials  and  Methods 

Clinical  Series 

We  satisfactorily  studied  151  of  155  consecutively  monitored 
major  trauma  patients  with  noninvasive  circulatory  monitoring 
beginning  shortly  after  their  admission  to  the  ED  and  continuing 
into  the  radiology  department,  the  OR,  through  the  postanesfoe- 
s-a  recovery  area,  and  to  the  ICU.  Four  patients  were  excluded 
because  of  insufficient  data  due  to  technical  equipment  failure  or 
communication  issues  with  personnel;  the  present  report  de¬ 
scribes  the  data  of  151  patients.  Table  1  lists  the  salient  clinical 
features  of  the  scries.  Patients  with  major  blunt  trauma  or 
penetrating  trauma  and  significant  risk  of  mortality  or  morbidity 
were  selected  for  monitoring  prior  to  possible  emergency  mr. 
g&ry.  The  criteria  for  resuscitation  were  empirically  determined 
by  the  previous  series  of  survivors1  values  and  by  the  best  possible 
mitlai  responses  of  the  cardiac  index  (Cl)  and  the  other  hemo¬ 
dynamic  variables. 'Z-K  Monitoring  was  continued  until  a  plateau 
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Table  1  —-Clinical  Feature** 


ATfunsurvivors 

-  VBriab1es  (n  =  4,1)  pValuct 

35  ±4  40  2  3  025 

Gender  ’ 

90  (3?)  41/48  (85)  0.937 

MS2mof  M81 ^  0937 

injury 

Gunshot  wound  42/61(69)  19/61(3.1)  0001 

Bhmt  tmuma  41/68(60)  27/68(40)  003 

S«b  wound  80/22(91)  2/22(9)  0001 

Bodily  injuity 

*?cad  24/47.(59)  17/41(41)  '•  0.759 

Chest  50/68(73)  18/65(26)  0.00J 

Abdomen  54/84  (64)  30/84(36)  0.001 

CCS  13.3  ±  0.3  9,4  2  0.7  0  00,7 

Injury  severity  score  27-8  2  4.7  30.5  ±  4.6  0.24 

♦Values  given  as  No.  (%)  or  mean  ±  SD  “  - 

tp  Values  for  differences  between  the  number  and  percentage!  of 
Survivors  vs  nonsurvivors. 


- ...  ....  pvKMure,  our  rne  adequacy  ol  initial  resusci- 

tafaons  may  have  been  limited  in  part  by  clinical  exigencies  at  tine 
time.  The  calculation  of  cumulative  excesses  or  deficits  and 
discriminant  analysis  were  performed  after  monitoring  wnscom- 
pleted.  The  institutional  review  board  approved  the  protocol. 

Mean  Arterial  BF 


Continuous  mean  arterial  BP  (MAP)  was  measured  noninva- 
stvely  (Dynamap  system;  Critieon;  Tampa,  FL)  or  was  calculated 
electronically  from  transducers  in  liue  with  intra-arterial  cathe¬ 
ters  when  the  latter  were  used. 


Cardiac  Output 

A  thoracic  bioelectric  impedance  device  (IQ  system;  Wnntagh 
Inc;  Bristol,  FA)  was  applied  shortly  after  the  arrival  of  the 
patient  in  the  ED.  Fairs  of  noninvasive,  disposable,  prewired 
hydrogen  electrodes  were  positioned  with  one  pair  placed  on 
each  side  of  the  base  of  the  neck  and  two  other  pairs  placed,  one 
on  each  side  of  the  chest  at  the  level  of  the  aypbisteriial  junction 
opposite  the  lateral  axillary  line.  Three  ECG  leads  were  placed 
across  the  precordium  and.  left  shoulder.^  A  1.00.KHz,4-m  A 
alternating  current  was  passed  through  foe  patient’s  thorax  by  foe 
Outer  pairs  of  electrodes,  and  the  voltage  was  sensed  bv  the  inner 
pairs  of  electrodes;  foe  voltage  sensed  by  the  inner  electrodes 
captured  foe  baseline  impedance,  foe  first  derivative  of  foe 
impedance  waveform,  and  foe  ECG.  The  ECG  and  bioimped¬ 
ance  signals  were  filtered  with  an  all-integer-cocfficieot  technol¬ 
ogy  to  decrease  computation  and  signal-processing  times  The 
signal  processing  algorithm  used  a  time-frequemy  distribution 
(modified  Wigner  distribution)  analysis  that  increased  sienal-to- 
noise  ratios.***!  The  data  were  automatically  acquired  and 
downloaded  to  a  floppy  disk.  When  indicated  by  clinical  criteria, 
FACs  were  inserted  into  the  patient:  in  the  OB  or  the  IGU  and 
CI  estimations  ware  made  at  least  hourly  in  unstable  patients  and 
cveiy  4  h  in  stable  patients.  The  optimal  goal  for  CI  in  various 
etiologic  diagnostic  groups  was  defined  by  survivors'  values '*-•■' 
axid  was  tested  in  subsequent  studies.1'1  -si 
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Limitations  of  the  impedance  method  include  faulty  electrode 
placement,  motion  artifacts,  restlessness,  shivering,  pulmonary 
edema,  pleural  effusion,  valvular  heart  disease,  dysrhythmias,  and 
electrical  leaks  from  other  instruments  using  the  same  circuit, 
These  are  usually  apparent  from  inspection  of  the  impedance 
waveform  and  by  the  following  previously  described  criteria; 
baseline  impedance  >  15  ohms  and  impedance  signal  >  0.3  ohm, 
which  usually  indicate  pulmonary  edema  due  to  cardiac  failure  or 
kte-stage  ARDS.11  These  limitations  were  excluded  during  the 
time  of  monitoring  in  the  present  study. 


False  Oximetry 

Arterial  oxygen  saturation  (Sao2)  was  assessed  continuously  by 
pulse  oximetry  (Nell cor;  Pleasanton,  CA)  as  a  reflection  of 
pulmonaiy  gas  exchange.  Values  were  observed  and  recorded  at 
the  time  of  the  Cl  measurements.  Appreciable  or  sudden 
changes  in  these  values  also  were  noted,  and  changes  to  <  94% 
were*  confirmed  by  SaOs  measurement  obtained  by  standard 
blood  gas  analysis.11-1* 

Transcutaneous  Oxygen  Tension 

Standard  transcutaneous  oxygen  tension  (tcFOa)  measure¬ 
ments  were  continuously  monitored  throughout  the  observation 
period.  This  technology  uses  the  same  Clark  polarograpbie 
oxygen  electrode  routinely  employed  in.  standard  blood  gas 
nioasureTncnts,^-*0  The  oxygen  tensions  were  measured  in  a 
representative  area  of  the  skvn  surface  heated  to  4*C  to  increase 
diffusion  of  oxygen  across  the  stratum  comeum  and  to  avoid 
vasoconstriction  in  the  local  area  of  the  skin  being  measured.*7 
Previous  studies  demonstrated  the  capacity  of  transcutaneous 
oxygen  tensions  to  reflect  tissue  oxygen  tension.11-12*2**®  tcPo2 
has  been  shown  to  reflect  the  delivery  of  oxygen  to  the  local  area 
of  skin,*  it  also  parallels  the  mixed  venous  oxygen  tension  except 
under  late  or  terminal  conditions  In  which  peripheral  shunting 
leads  to  .high  mixed  venous  hemoglobin  saturation  values*'' While 
oxygen  tension  of  a  segment  of  the  skin  does  not  reflect  .the  state 
of  oxygenation  of  all  tissues  and  organs,  the  skin  has  the 
advantage  of  being  the  most  sensitive  early  warning  tissue  of  the 
adrenomedidlaiy  stress  response;  vasoconstriction  of  the  skin  is 
an  early  stress  rehouse  to  hypovolemia  and  other  shock  *yn- 
dromes.11’13-**  tcFOz  Values  were  indexed  to  the  fraction  of 
inspired  oxygen  (Fio?i)  concentration  to  give  a  taPc^FlOa  ratio 
because  0f  parked  tcPo2  changes  produced  by  changes  in  the 
level  of  inspired  oxygen.  The  thermal  environment  was  maintained 
at  reasonably  constant  levels,  and  marked  changes  in  room  temper¬ 
ature  tom  drafts  or  open  windows  were  avoided  to  maintain  the 
accuracy  of  the  transcutaneous  methods.  In  addition,  the  electrode 
must  be  moved  to  a  nearby  thoracic  or  shoulder  site  every  4  h  and 
recalibrated  to  avoid  first-degree  skin  bums- 

Level  of  Consciousness 

At  the  time  of  the  patient's  admission  to  the  ED,  the  clinical 
team  evaluated  and  recorded  the  degree  of  unconsciousness  by 
the  Glasgow  coma  scale  (GCS).  which  uses  aye  movement,  verbal 
responses,  and  motor  responses  to  verbal  and  painful  stimuli.  The 
clinical  service  also  noted  changes  in  the  GCS  throughout  the 
patient’s  hospital  course. 

Estimated  Blood  Loss  ett  the  Time  of  Surgery 

Blood  loss  was  estimated  hy  die  surgeon  and  anesthesiologist 
intraoperatively  in  a  routine  manner  by  counting  lap  tapes  and 
sponges  and  by  measuring  the  contents  of  suction  botdes. 


Method  for  Calculating  the  Total  Cumulative  Excess  or  Deficit 
of  Each  Monitored  Variable 

The  patterns  of  each  patient  were  examined  for  motion  artifact, 
noise,  effect*;  of  fluid  and  vasopressor  therapy,  manipulation  of 
tubing,  and  other  extraneous  factors.  The  total  overall  deficit  or 
excess  of  each  noninvasively  monitored  variable  was  evaluated  by 
comparing  its  normal  or  optimal  value  with  its  temporal  pattern 
during  the  oblation  period.  This  was  done  by  mathematically 
integrating  over  time  the  area  between  the  continuous  display  of 
each  fluctuating  variable  and  either  die  normal  values  for  BP, 
Sao*  and  tcPo/ftoa  or  the  optimal  goal,  as  defined  by  the  Cl 
values  of  Survivors  during  the  first.  24  h  after  hospital  admis¬ 
sion.11-*1 

The  net  cumulative  deficits  or  excesses  were  calculated  for 
each  individual  patient  and  for  both  survivor  and  non  survivor 
groups  as  time-integrated  areas  between  the  curve  produced  by 
continuously  monitored  variables  and  their  normal  or  reference 
values.  For  example,  given  a  normal  MAP  of  85  mm  Hg,  in  a 
patient  whose  MAP  averaged  60  mm  Hg  for  2  h  before  resusci¬ 
tation*  tbe  calculated  deficit  is  -50  ram  Hg/h  ([85-60]  *  2). 

Flow  calculations,  measured  as  volume  per  unit  of  time,  are  in 
liters  per  minute  per  square  meter.  When  multiplied  by  the 
monitored  time  in  minutes,  this  gives,  as  units,  liters  per  square 
meter  for  01  or  liters  for  cardiac  output  The  units  for  MAP,  Sac* 
and  tcPojj/FiOa  arc  millimeters  of  mcrcuiy  per  hour,  percent  per 
bout,  and  millimeters  of  mommy  per  hour,  respectively. 

When  the  mean  MAP  deficits  were  calculated  using  all  values, 
a  large  number  of  normal  high  values  obscured  the  deficits;  the 
patients  with  cardiac  arrest  and  zero  MAP,  for  example,  showed 
no  net  MAP  deficit,  because  the  many  normal  and  high  values 
overshadowed  the  later  short  but  lethal  hypotensive  episode.  For 
MAP,  therefore,  we  calculated  cumulative  deficits  from  de¬ 
creases  below  the  normal  range. 


Statistical  Analysis 

The  survivors 1  andnonsurviyors’  deficits  of  MAP,  Cl,  Sao2,  and 
tcPOg/Flo*  were  calculated  for  the  periods  of  monitoring-  Each 
of  die  categoric  variables  was  tested  for  the  difference  in 
distributions  between  the  two  outcome  groups,  those  who  sur¬ 
vived  and  those  who  died  during  the  current  hospitalization, 
using  the  x2  test  or  two-tailed  Fisher's  Exact  Test.  The  t  test  with 
Bonferronj  correction  was  applied  to  each  of  the  continuous 
variables  to  compare  the  means  of  the  two  outcome  groups. 
Variables  ^considered  for  discriminant  analysis  were  Cl,  GCS, 
SaoBI  tePo2/.FTo2,  MAP,  heart  rate,  Pao2,  hematocrit,  transcuta¬ 
neous  C02  tension  (PtcC02),  injury  severity  score^age,  and 
gender.  The  first  four  met  the  criteria  (p  <  0.20).  *  * 

The  variables  that  were  significant  at  the  p  <  0.2  level  by  the 
aforementioned  x2  tests  or  the  t  tests  were  fed  into  a  stepwise 
discriminant  analysis  (PR.OC  STEPdisk)  to  identify  the  variables 
that  collectively  contribute  to  differentiate  th«  two  outcome 
groups.  Thus,  the  variables  selected  then  were  entered  into  a 
model  in  PROC  DISCRJ.M  to  derive  the  discriminant  function 
by  generalized  squared  distance,  taking  into  account  the  prior 
probabilities  of  the  groups.  This  procedure  evaluated  the  dis¬ 
criminant  function  by  calculating  the  error  rate  estimates  or  the 
probabilities  of  miscJassification. 

Cross-validation  of  the  results  was  performed  by  the  jackknife 
method.  The  data  were  split  into  two  independent  samples  by 
taking  the  data  of  every  other  patient.  One  group  was  used  for 
calibration  to  generate  another  series  of  classification,  functions, 
and  the  remaining  group  was  used  to  calculate  results  based  on 
the  new  classification  functions.  The  statistical  analyses  were 
performed  with  a  computer  program  (SAS  for  Windows,  Release 
6.12;  SAS  Institute;  Gary,  NC). 
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—  -  Survive rs 
Nonsurvivors 

Ftcube  1.  The  temporal  patterns  of  survivors  (solid  line)  and 
nonsumvoTs  (dashed  line?  for  MAP,  Cl,  $aC2  (Sap02),  and 
tcPo2  (Ptc08)  indeed  to  the  tcPOj/Fio2  ratio.  AlJ  values  axe 
keyed  to  the  time  of  admission  to  the  ED.  Dote  represent  mean 
vaives,  and  vertical  lines  represent  SEM.  Crow-hatched  arpos 
indicate  the  normal  range  for  MAP,  Sao*  and  tcPo*/FTOE  ratio 
and  optimal  goals  for  Cl  Note  that  the  Cl,  MAPrSao!  and 
tcPo/Fio*  values  of  survivors  were  generally  higher  chan  those 
of  the  noasuTvivors. 


Results 

Noninvasive  Monitoring  From  the  Time  of 
Admission 

_  Thf  use  of  noninvasive  monitoring  systems  was 
foiina  to  be  feasible  in  patients  experiencing  short¬ 
term  emergency  conditions  for  the  early  description 
of  temporal  hemodynamic  patterns  and  to  provide 
quantitative  calculations  of  the  total  amount  of  def¬ 


icit  or  excess  accumulated  by  each  monitored,  vari¬ 
able.  There  were  103  survivors  and  48  nonsurvivors 
(mortality  rate,  32%),  131  patients  were  men,  20 
patients  were  women,  and  the  average  (±  SEM)  age 
was  —  1.4  years.  Of  61  patients  with  gunshot 
wounds,  42  survived  and  19  died  (mortality  rate, 
31%).  Of  68  patients  with  blunt  trauma.  41  survived 
and  27  died  (mortality  rate,  40%).  Of  22  patients  who 
sustained  stab  wounds,  20  survived  and  2  died 
(mortality  rate,  9%).  Of  41  patients  who  had  head 
injuries,  24  survived  and  17  died  (mortality  rate. 
41%).  Of  68  patients  who  sustained  chest  injuries,  50 
survived  and  18  died  (mortality  rate,  26%).  Of  84 
patients  who  sustained  abdominal  injuries-,  54  sur¬ 
vived  and  30  died  (mortality  rate,  36%).  Sixty-eight 
patients  had  injuries  involving  more  than  one  bodily 
area.  The  injury  severity  'score  (±  SEM)  was 
21.8  ±  4.7  for  survivors  and  30.5  ±  4.6  for  nonsur¬ 
vivors  (p  =  0.24), 

Monitoring  was  performed  for  7.9  ±  2.6  h  during 
the  initial  resuscitation  (survivors,  7.8  h;  nonsurvi¬ 
vors,  8,3  h).  Subsequently,  survivors  were  monitored 
intermittently  to  15.6  ±  7.1  h  after  hospital  admis¬ 
sion,  and  nonsurvivors  were  monitored  to  1.8.7  ±  8.4 
h  after  hospital  admission. 

.  The  data  of  emergency  patients  from  the  time  of 
their  ED  admission  are  shown  in  Figure  1.  The 
correlation  between  simultaneous  thermodilution 
and  bioimpedance  cardiac  output  measurements  in 
the  present  series  was  r  =  0.91  and  r2  -  0.83,  and 
bias  and  precision  were  -0.30  ±  1.10  L/min/m2. 
Table  2  lists  the  mean  ±  SEM  of  Cl,  MAP,  Sao2, 
and  tcPo^/FlOa  for  survivors  and  non$urvj.vors  aver- 
aged,  throughout  the  observation  period.  The  Cl, 
SaOa,  and  tcPoa/FlOg  values  of  patients  who  survived 
were  significantly  greater  than  for  those  who  .died. 
MAP  values  of  survivors  tended  to  be  higher  than 
those  for  nonsurvivors  (p  =  0.066)  (Table  2). 

The  body  temperatures  of  survivors  and  nonsurvi- 
vors  at  hospital  admission  averaged  36.7  ±  Q,0°C  to 
7°C  and  36.3  ±  0.0°C  to  9°C,  respectively.  'We  took 
aggressive  precautions  to  correct  hypothermia  when 
it  occurred,  especially  in  the  OR  where  conditions 
were  more  controllable. 


Table  2 — Noninvasive  Hemodynamic 

Value*  for  Survivor*  and  Nonxurvivon* 

NoimaJ  or 

Variable  Optimal  Value 

Survivors 
(n  -  103) 

Nonsurvivors 
(n  -  48) 

p  Vo.Iuel 

Cl,  4  0 

Map,  tom  Kg  05 

Sno2,  % 

tcPOg/Flo^,  mm  Kg  2O0 

4.14  ±  0.02 

88  ±  0.37 

99  ±  0.05 

206  £  2,0 

3.87  ±  0.03 

80  £  0.69 

96  ±  0.26 

93  ±  2.6 

<  0.001 
0.068 

<  0.001 

c  0  nni 

-values  given  aa  meat*  S  SEM,  unless  otherwise  indicated,  - - - - - __ 

tp  Values  for  differences  between  survivors'  and  nonsuivivms'  values. 
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Table  3 — Time  to  Reach  Coal  in  Patients  Who  Attained  Coal * 


Venable 

Survivors 

Nonsurvivors 

1 

Hours 

No J%  Goals' 

Not  Beached  -* 

I 

Hours 

No,/%  Goals 
Not  Reached 

Cl 

1.99  ±  5,10 

6/6 

2.64  ±  3.69 

13/28 

MAJ> 

1,34  ±  2,80 

0/0 

2.03  ±  3.69 

16/34 

Sjo^j  - 

0.48  X  1.66 

o/o 

1.86  ±  4.44 

9/19 

tdPOa/FlOjj 

3.14  ±  4.61 

_ * - L - - - 

7/7 

4-48  ±  4.56 

39/83 

♦VaJues  given  as  mean  ±  $D,  unless  otherwise  indicated. 


The  mean  (±  SD)  estimated  blood  loss,  which  re¬ 
flects  preoperative  and  intraoperative  hemorrhaging, 
measured  2,970  ±3,856  mL  in  survivors  and 
6,263  ±  5,540  mL  in  the  nonsurvivors  at  the  end  of 
surgery.  In  the  present  series,  there  were  22  patients 
wb-0  had  massive  blood  loss  (ie,  >  5,000  mL).  Vigorous 
attempts  were  made  to  replace  these  losses  at  the  time 
of  surgery  and  in  the  immediate  postoperative  period. 

Temporal  Circulatory  Patterns  in  Survivors  and 
Nonsurvivors 

Figure  1  shows  the  temporal  patterns  of  noninva- 
sive  circulatory  variables  of  the  survivors  and  non¬ 
survivors  beginning  with  the  initial  measurements 
after  admission  to  the  ED.  Cl  values  were  initially 
higher  in  the  survivors.  The  Sao2  values  of  nonsur- 
viv'ors  were  significantly  lower  than  the  those  of 
survivors,  but  these  differences  were  not  clinically 
important;  when  Sao2  reductions  occurred,  they 
wesre  rapidly  corrected  by  intubation,  mechanical 
ventilation,  or  increased  Fio2.  The  values  for  the 
tcPo/Flo2  ratios  of  nonsurvivors  were  markedly 
lower  than  those  of  survivors  and  were  lower  than 
normal  throughout  the  observation  period.  Table  3 
lists  the  time  taken  to  achieve  goals  of  therapy  for 
ea<h  variable  that  reached  the  desired  end  point  as 
well  as  the  number  and  percentage  of  those  who  did 
not  reach  the  goals.  The  deaths  of  nonsurvivors 
occurred  an  average  of  8.7  ±  2.8  days  after  hospital 
admission.  However,  there  was  a  bimodal  distribu¬ 


tion  with  17  deaths  in  the  first  8  h  and  14  deaths 
occurring  ^  10  days  after  hospital  admission. 

Net  Cumulative  Amount  of  Deficit  or  Excess  in 
Monitored  Variables 

Table  4  shows  the  net  cumulative  deficit  or  excess 
of  monitored  variables  used  to  evaluate  cardiac, 
pulmonary,  and  tissue  perfusion  functions.  Figure  2 
is  an  illustrative  example  of  a  survivor  whose  Cl  and 
tissue  perfusion  deficiencies  were  corrected  at  19 
and  23  h  postadmission,  respectively.  Figure  3  shows 
the  data  of  a  patient  whose  Cl  and  tissue  perfusion 
deficiencies  persisted  for  >  24  h.  He  developed 
lethal  ARDS. 

Outcome  Prediction 

There  were  significantly  greater  calculated  deficits 
of  Cl,  pulse  oximetry,  and  transcutaneous  02  in 
nonsurvivors  than  in  survivors  during  the  period  of 
monitoring  (Fig  4  and  Table  4).  These  three  vari¬ 
ables  and  the  GCS,  having  moderate  levels  of  signif¬ 
icance  with  outcome,  were  selected  for  the  stepwise 
discriminant  analysis  (PROC  STEPdisk).  Based  on 
the  classification  ■  function  generated  for  each  of 
these  four  variables  in  PROC  DISCRIM,  the  dis¬ 
criminant  function,  Z,  was  derived:  ^  ^ 

Z  =  O.OOlla-b  0.3300b  +  0.0656c  +  0.0423d 
where  a  represents  cumulative  tcPp2/Flo2  values,  b 


Table  4 — Mean  Net  Cumulative  Deficits  or  Excesses  of  Monitored  Values  of  Survivors  and  Nonsurvivors 

Throughout  the  Period  of  Observation  _ _ 


Survivors 

Nonsurvivors 

Variable 

r . . . 

Mean 

SEM 

Mean 

SEM1 

p  Value* 

cr.  Urn2 

M-AJP,  mm  Hg/h 

Sao2,  %/h 

tc^Oj/FiOsi,  mxn  H^!h 

+  81 
-  10 
-  1 
+  313 

52 

13 

0.3 

87 

iiii 

138 

24 

2.6 

175 

<  0.007 
0.078 

<  0.006 
<  0.001 

*j>  Values  are  for  differences  between  survivors’  and  nonsurvivors’  values. 
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represents  the  GCS,  c  represents  cumulative  Sao2 
values,  and  d  represents  cumulative  CT  values.  Table 
5  summarizes  the  relative  influence  of  each  variable 
with  respect  to  outcome.  Ninety-five  percent  of  the 
survivors  and  62%  of  the  nonsurvivors  were  correctly 
classified  in  the  first  24  h  postadmission  (Table  6),  Of 
151  patients,  23  (15.2%)  were  misclassified.  Five  of 
the  35  patients  predicted  to  die  in  the  first  24  h 
subsequently  improved  and  lived. 

Results  of  Cross-Validation 

Cross-validation  of  the  discriminant  analysis  by  the 
jackknife  method  demonstrated  results  that  were 
similar  to  the  initial  calculation  for  the  series  as  a 
whole.  The  results  of  the  calibration  data  set 


(N  =  75)  are  shown  in  Table  7,  and  the  result?  from 
thevalidation  data  set  (N  =  76)  are  shown  in  Table 
8.  The  cross-validated  discriminant  analysis  was 

Z  +  0.0018a  +  0.3138b  +  0.0786c  +  0.0022d 

where  a ,  b,  c,  and  d  are  defined  as  above.  The 
classification  of  die  survivors  was  Z  >  1.91,  Miscal¬ 
culations  occurred  in  12  of  76  patients  (16%).  This 
was  considered  to  be  in  satisfactory  agreement  with, 
foe  initial  calculation  and  suggests  consistency  of  the 
data  by  this  analysis. 

Discussion 

The  limitations  of  noninvasive  bioimpedance  car¬ 
diac  output  monitoring  include  motion  artifacts, 
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FIGURE  3.  Data  from  a  26-yRar-old  man  who  sustained  multiple  stab  wounds  of  the  abdomen  with 
lacerations  of  the  stomach,  duodenum,  and  superior  mesenteric  vein.  He  bod  marked  reduction  of  cardiac 


to  dobutaminc  and  dopamine  inhisions  for  an  estimated  18?0b0-mL  blood  loss.  Hie  patient  died  < 
and  multiple  organ  failure.  See  tbe  legend  of  Figure  1  for  abbreviations  not  used  in  the  text  . 


arrhythmias,  pulmonary  edema,  pleural  effusions, 
and  expansion  of  interstitial  fluid  from  massive  crys¬ 
talloid  infusions.  The  advantages  of  this  noninvasive 
monitoring  system  include  technical,  convenience 
and  the  continuous  display  of  data  allowing  the 
calculation  of  the  amount  of  deficit  or  excess  of  each 


variable,  from  the  time-integrated  area  under  the 
curve.  The  area  under  the  curve  provides  an  arith¬ 
metic  solution,  to  replace  the  subjective  evaluation  of 
irregular  curves  and  provides  estimates  of  cardiac, 
pulmonary,  and  tissue  perfusion  functions. 


Table  5 — Steptrise  Discriminant  Analyst** 


Step  Entered 

Partial*2 

Prob  >  P 

Cumulative 

Cumulative  tcPOa/F.TUa 

0.210 

0.0001 

0.2000 

%  GCS 

0,188 

0.0001 

0.3581 

3.  Cumulative:  Sao2 

0.053 

0.0047 

0,3021 

4.  Cumulative  Cl 

0.031 

0.0336 

0.4107 

^Classification  of  survivors,  Z  >  2.3(3;  where  Z  =  0.0011  (cumulative 
tcPo^/FtC^)  +  <1.3300  (GCS)  +  0.0656  (eirniulartw  Sao2)  +  0.0423 
(cumulative  Cl). 

fPfllaTs  trace/(No.  of  groups  —  1). 


Table  6 — Classification  Summary  for  the  Serins 
(a  =  151)* 


Predicted 
to  Die 


Predicted 
to  Live 


Outcome 

i - 

No. 

- — I 

Rcrw% 

l - 

No, 

i 

B.ow% 

No. 

Cot^' 

Died 

30 

62.5 

IS 

37.5 

48 

3.1.8 

Lived 

5 

4.9 

98 

95,1 

103 

68,2 

Total,  % 

35 

23.2 

116 

76,8 

151 

100.0 

♦Misclassification;  23/151  (1  $&%).  Row% 

=  the  percentage  oF  patients 

In  that  row;  CoI%  =  the  percentage  of  patients  in  that  cotumn. 
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Table  7 — Classification  From  the  Calibration  Dataset 
(n  =  75)* 


Table  8— Classification  From  the  Validation  Dataset 
(n  =  76)* 


Predicts?  Predicted 


Actual 

to  Die 

to  Live 

Totul 

Outcome 

No, 

.  Row% 

i 

No. 

Buw% 

i - 

No. 

Co!% 

Died 

15 

(32.5 

9 

37.5 

24 

32.0 

Lived 

3 

5.9 

4tS 

94.1 

51 

(38.0 

Total,  % 

18 

24,0 

57 

76.0 

75 

100.0 

'Misckssificafcion:  12/75  (16.0%),  See  Table  6  for  abbreviations  not 
used  in  text. 


The  net  cumulative  deficits  of  flow  and  tissue 
perfusion  measured  during  the  initial  resuscitation 
period  were  greater  in  uonsurvivors  than  survivors; 
these  differences  were  correlated  with  outcome.  For 
example,  during  the  monitoring  period,  the  Cl  val¬ 
ues  of  survivors  averaged  8.1  L/m2  more  than  the 
optimal  4.0  L/min/m2,  which  was  determined  empir¬ 
ically  from  the  plateau  of  high  values  of  survivors 
within  the  first  24  h  of  hospital  admission.11  21  This 
was  equivalent  to  140  L  of  cardiac  output  per  patient 
over  die  monitored  period.  During  the  monitoring 
period  of  those  who  died,  the  Cl  averaged  232  L/m5 
less  tli an  optimal,  and  the  cardiac  output  averaged 
402  L  per  patient  less  than  optimal.  The  difference 
between  survivors  and  nonsurvjvors  was  542  L.  We 
used  4.0  TVmiu/ra“  as  the  therapeutic  goal  because 
this  was  the  mean  value  for  the  first  24-b  period,  on 


Predicted  Predicted 

Actual  t0  toLiv<i  Total 

Outcome  No.  Row%'  No,  Row%'  Wo.  Cd% 

Died  15  62.5  9  37.5  24  3H5 

Lived  5  9.6  47  90.4  52  fifl.4 

Total,  %  20  26.3  56  73.7  76  100.0 

* M isclnssl fico tioil;  14/76  (18%).  See  Table.  6  for  abbreviations  not 
used  in  text. 


which  this  study  was  focused.  This  goal  admittedly  is 
arbitrary  and  points  to  the  need  for  additional  re¬ 
search  in  this  area. 

The  high  early  Cl  values  in  survivors  suggest  that 
there  may  have  been  less  hypovolemia  and/or  better 
physiologic  compensations.  This  concept  is  rein¬ 
forced  by  the  greater  tcPOj/FiOz  net  cumulative 
excesses,  which  suggest  better  tissue  perfusion/oxy- 
genation  for  survivors  in  the  initial  stages.  These 
preliminary  studies  need  to  be  evaluated  indepen¬ 
dently  in  larger  series  with  different  types  of  acute 
illnesses  and  emergency  conditions.  Furthermore, 
additional  studies  are  needed  to  evaluate  the  effects 
of  specific  trunk  and  extremity  traumas,  head  inju¬ 
ries,  pelvic  and  long  bone  fractures,  prior  organ 
dysfunctions,  and  other  comorbid  states  on  the 
validity  of  this  oarly  predictive  model. 


CI,L/m2 


FIGUSE  4.  Cumulative  excesses  and  deficits  in  survivors  and  nonsurvivors  for  CT  and  tcPOo/FlO, 
calculated  tor  the  monitored  period.  See  the  legend  of  Figure  1  ior  abbreviations  not  used  in  the  text 
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•  underlying  this  approach  is  that 

circulatory  deficiencies  that  ultimately  lead  to  shock 
organ  failure,  and  death  may  be  identified  eatt; 
nomnvasive  monitoring  even  in  the  extenuating 
cumstances  of  severely  traumatized  emergency  pa- 

dia^w  alar2e  mner  city  public  hospital  Earlier 
diagnosis  of  a  circulatory  deficiency  allows  therapy  to 
be  initiated  sooner  in  the  hope  that  earlier  therapy 

crucid1Pr°Ve  °UtCOme  in  ernergencies  where  time  is 
More  importantly,  noninvasive  monitoring,  which 

sensitive,1^*  allows  estimates  of  the  amount  of 
deficits  calculated  from  the  difference  in  the  areas 
between  normal  values  or  survivor  values  and  the 
continuously  monitored  variables.  Multiple  noninva¬ 
sive  hemodynamic  monitoring  systems  provide  sim- 

inm’f’1"’3?"  *°  ?“  °f  FAC,  except  for 
p^naiy  ^ty  ocdMjfon  pressures.  Discrtofoant 
analysis  of  these  data  provides  a  mathematical  basis 
for  outcome  prediction.  Future  prospective  clinical 

SZ, a!  0ther  rafdtutions  nceded  to  validate  the 
present  approach. 

3,80  Provide5  “  approach 
&  TGi  t0  j6Ve  °P  ™  organized  coherent 
therapeutic  plan  based  on  physiologic  criteria  for  the 
^mergency  patient  as  he/she  proceeds  from  the  ED 
o  the  OR,  the  radiology  department,  and  the  ICU 
Linear  discriminant  function  predicted  outcome  cor¬ 
rectly  m  95%  of  the  survivors  and  in  62%  of  the 

in  after  hosPitai  aduus- 

sion-  This  was  probably  as  much  as  should  be 

expected  for  nonsurvivors  since  many  patients  devel- 
oped  lethal  complications  unrelated  to  their  injuries 
late  m  their  hospital  course, 

Smce  the  essence  of  tissue  perfusion  is  an  ade¬ 
quate  supply  of  oxygenated  blood  to  the  tissues 
perfusion  is  Inferred  from  the  direct  measurement  of 

ithrSgenatiQn  !?"«  **  Clark  pokrographic 
method  for  oxygen  tension.^  Although  the  skin  is 

°^PTentative  ofa®  Hssues> is  *6  largest  organ 
and  the  first  organ  to  be  affected  by  the  achenomed- 
ufiary  stress  response.  tePo2  provides  early  warning 
in  acutely  rll  emergency  patients^-  it  tracks  oxygen 
uptake  in  acute  clinical  shock  episodesU  and  JS 

°LexPerimental  hemorrhage 

oaSac  '  U  Ca?aC  ^  resP^tory  Mure, 

^  I  f  ^^'^ary  resuscitation  in 

A.  Shown  in  the 

%S m*mre of  ttsSM p-*— 

common  term,  such  as  the  combination  of  Cl  md 


3^! were  avoid®d  to  minimize  statistical 
problems  of  discriminant  analysis.  This  does  not 
mean  that  the  more  conventional  variables  like 
tachycardia,  hypotension,  acidosis,  skin  color,  lactate 
fevek  mental  status  etc,  are  not  usefoj  at  *n£ 
hen  they  occur.  Obviously,  when  they  are  abnor¬ 
mal,  they  are  extremely  useful  and  important.  How- 
ever  the  criteria  of  the  present  study  focused  on 
earfy  nomnvasive  hemodynamic  variables  in  the  im¬ 
mediate  postadmission  period  that  most  consistently 
separated  survivors  and  nonsuivivors.  y 

The  concept  that  hypovolemia  is  an  early  primary 
problem  that  plays  an  important  role  in  low; flaw  and 
poor  tissue  perfusion  states  is  supported  by  the 

2STEi  (  idireCa  #""■*»  of^sive  hemor- 
(2>  estimated  blood  loss  of  hemoperitoneum 
and  hemothorax  at  the  time  of  surgeiy  in  patients 
who  underwent  surgical  exploration;  and  (3)  prior 
JSJ"  f  Ae  ?torature  that  documented  Wood 
vohime  deficits  m  posttraumatic  and  postoperative 

and  ched^°  Subse<IuentI>r  developed  organ  failures 
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